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Preface 


This  paper  culminates  a  nearly  five  year  effort  in  pursuit  of  a  PhD  degree  on 
a  part  time  basis.  The  last  two  years  were  spent  on  this  particular  research.  The 
primary  aim  of  this  study  was  to  propose  a  new  high  resolution  spectroscopic 
technique  using  a  nonresonant  field  to  effect  a  level  crossing.  The  theoretical 
description  of  the  Townes— Merritt  effect,  quantum  beats,  and  level  crossings  are 
covered,  as  well  as  a  brief  look  at  some  of  the  important  experimental  parameters. 

An  actual  experiment  is  suggested  as  a  logical  follow-on  effort.  A  secondary 
objective  was  to  gain  a  more  detailed  education  in  quantum  electrodynamics  (QED) 
and  nonlinear  effects  such  as  level  crossings  and  quantum  beats.  Several  useful 
techniques  in  atomic  physics  are  covered  in  this  paper,  which  hopefully,  the  reader 
may  find  useful. 

There  are  several  people  to  whom  I  owe  many  thanks  in  the  completion  of 
this  program.  First,  I  wish  to  thank  Dr.  Richard  J.  Cook  who  has  both  inspired  me 
in  this  field  of  study  and  has  ultimately  made  it  possible  for  me  to  pursue  and 
complete  this  program.  It  has  been  an  absolute  pleasure  to  work  with  and  learn 
from  him.  I  also  owe  deep  appreciation  to  Dr.  Ronald  L.  Bagley  for  his  willingness 
to  support  my  desire  to  take  on  such  a  program  and  constant  urging  to  complete  it. 

Finally,  and  most  of  all,  I  wish  to  thank  my  family,  Kristen,  Lara,  and  particularly 
my  wife  Barbara,  for  their  continued  support  and  patience  in  this  often  trying  task. 

It  is  to  them  that  I  truly  owe  this  opportunity  to  attain  a  lifelong  goal.  _ 
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Abstract 

A  high  resolution  spectroscopic  technique  using  a  nonresonant  RF  field  in  a 
level  crossing  experiment  is  proposed.  Absorption  sidebands,  originally  observed  by 
Townes  and  Merritt  (1947),  are  analyzed  using  a  quantum  electrodynamic  (QED) 
approach.  A  model  is  developed  using  the  Heisenberg  operator  formalism  and  the 
adiabatic  approximation  to  describe  the  atomic  dynamics.  Results  agree  with  a 
dressed  states  approach  and  the  original  experiment.  The  model  is  extended  to  a 
three  level  atom  to  describe  level  crossings  and  quantum  beats  for  linear  and 
quadratic  Stark  shifted  atoms.  Expressions  are  developed  to  define  level  crossing 
conditions  and  estimate  critical  experimental  parameters. 
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THEORETICAL  STUDY  OF  THE  TOWNES-MERRITT  EFFECT 
IN  LEVEL  CROSSING  EXPERIMENTS 

I.  Introduction 

The  aim  of  this  study  is  to  develop  the  theory  to  describe  the  effect  of  an 
external  nonresonant  electromagnetic  field  on  level  crossing  experiments.  The 
primary  motivation  rests  in  the  possibility  of  using  these  effects  to  obtain  greater 
accuracy  in  high  resolution  spectroscopy.  The  basis  for  such  a  study  began  with  the 
observation,  by  Townes  and  Merritt  in  1947,  of  sidebands  in  the  absorption  spectra 
of  the  OCS  J  =  1— >2  transition  when  the  sample  was  placed  in  a  nonresonant 
oscillating  electric  field  [1].  These  sidebands,  shown  in  Figure  1,  occur  at  intervals 
equal  to  twice  the  modulating  field  frequency.  The  original  semiclassical  analysis 
uses  time-dependent  perturbation  theory  to  describe  the  frequency  dependence  and 
relative  amplitudes  of  the  sidebands  [2]  with  good  agreement  to  the  experiment. 

The  proposed  application  of  this  "Townes-Merritt  Effect”  involves  the 
determination  of  energy  differences  between  nearby  atomic/molecular  levels  by 
causing  respective  sidebands  to  overlap  or  cross. 


Frequency  (MHz) 

Figure  1.  Absorption  sidebands  in  a  combined  static  and  RF  Field  [2(>] . 
Theoretical  pattern  is  below  experimental  result. 
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Level  crossings  and  a  closely  related  phenomenon,  quantum  beats,  can  occur 
when  there  exist  two  nearby  energy  levels  which  can  each  transition  to  a  lower 
state.  Both  result  from  interference  between  the  fluorescence  patterns  of  the  two 
independent  transitions  when  they  are  coherently  excited  with  a  short  pulse. 
Quantum  beats  manifest  themselves  in  an  oscillation  of  the  total  fluorescence  at  a 
frequency  proportional  to  the  energy  difference  of  the  excited  states.  Level  crossings 
result  in  a  rotation  of  the  overall  fluorescence  pattern  which  may  hp  detected  by  a 
change  in  intensity  or  polarization  in  a  given  direction.  These  will  be  discussed  in 
more  detail  later;  but  it  is  worth  mentioning  that  an  adequate  description  of  the 
phenomenon  requires  a  quantum  electrodynamic  (QED)  model  of  the  interaction 
[3,4],  which  will  become  the  basis  of  this  study. 

Level  crossing  experiments  are  an  established  method  used  to  measure  the 
fine  or  hyperfine  structure  of  atomic  and  molecular  systems  [5].  However  in  these 
cases,  static  fields  (electric  or  magnetic)  are  used  to  cause  the  level  crossing. 
Hereafter,  atomic  and  molecular  systems  will  be  referred  to  as  simply  atomic. 
Individual  levels  split  as  a  result  of  the  Stark  or  Zeeman  effect  and  orientation 
quantization  of  the  atomic  angular  momentum.  Then,  depending  on  the  field 
strength,  a  level  crossing  may  occur.  The  indication  of  a  level  crossing  is  a  change 
in  the  fluorescence  pattern  for  the  transition  from  the  crossing  levels  to  a  lower 
level.  The  original  energy  level  spacing  can  be  determined  from  the  field  strength  at 
which  the  crossing  occurs.  In  the  use  of  static  fields,  key  problems  present  are: 
strong  field  strengths  are  generally  required,  it  is  difficult  to  generate  static  fields 
with  uniform  strength  over  the  sample;  and  further  measurement  of  the  field 
strength  can  be  a  significant  source  of  error  [6j. 

Use  of  an  external  nonresonant  oscillating  field,  for  example  in  the  radio 
frequency  (RF)  or  microwave  region,  would  appear  to  alleviate  the  problems 
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I2/H-1  > 


Frequency  — * 

Figure  2.  Energy  Level  Diagram  for  a  combined  atom— RF  field  system. 

associated  with  static  fields.  The  primary  difference  in  the  case  of  an  RF  field  is 
that  the  level  splitting  is  dependent  on  the  external  field  frequency  instead  of  field 
intensity.  Figure  2  shows  a  typical  energy  level  diagram  of  an  atom  in  an  RF  field 
where  the  levels  represent  the  combined  atom— field  system.  Note  the  vertical  scale 
is  relative  to  the  energy  of  the  n  photon  state.  The  atom  with  no  field  applied  has 
two  closely  spaced  upper  states,  1 1  >  and  |2  ),  each  of  which  has  an  allowed 
transition  to  the  ground  state  jO  ).  As  the  field  is  applied  with  increasing 
frequency,  the  energy  levels  split  and  eventually  cross.  The  notation  used  to 
describe  specific  combined  states  indicates  the  original  atomic  state  and  the  photon 
number  state  n,  for  example  j  l,n).  It  is  worth  noting  that  for  each  atomic  state 
there  are  an  infinite  number  of  combined  states  due  to  the  infinite  number  of 
possible  photon  number  states;  however,  population  of  these  states  depends  on  the 
field  intensity.  Hence,  the  strength  of  givpn  transitions  and  observed  level  crossings 
would  depend  on  the  field  strength. 
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Figure  3.  Experimental  Setup 


The  potential  advantages  of  the  RF  field  approach  result  from  the  energy 
level  dependence  on  field  frequency.  The  level  crossing  points  would  depend 
primarily  on  field  frequency  which  can  be  more  accurately  measured  than  field 
strength.  This,  in  turn,  would  result  in  a  more  accurate  determination  of  the 
original  energy  level  separation.  The  major  impact  of  field  strength  would  be  in  the 
intensity  of  the  transitions.  Secondly,  unlike  the  static  case,  the  field  uniformity  is 
not  a  critical  concern  thus  eliminating  the  principal  source  of  error. 

The  proposed  technique  would  involve  an  experimental  setup  similar  to  the 
schematic  shown  in  Figure  3.  An  appropriate  sample  placed  in  a  variable  RF  field 
is  illuminated  by  a  pulsed  resonant  probe  beam  to  coherently  excite  the  transitions. 
A  detector  placed  off  axis  measures  the  fluorescence  of  the  sample  in  order  to 
observe  quantum  beats  or  to  indicate  a  level  crossing.  The  frequencies  at  which 
level  crossings  occur  determine  the  original  energy  level  spacing. 

The  main  objective  of  this  research  is  to  develop  the  theory  describing  the 
resonant  fluorescence  of  an  atomic  or  molecular  sample  as  a  function  of  the 
frequency  of  an  external  nonresonant  field.  Of  particular  interest  is  a  situation  with 
closely  spaced  excited  states  which,  due  to  the  Townes-Merritt  effect,  split  and 


eventually  result  in  level  crossing.  The  original  spacing  can  then  be  determined  by 
the  frequency  at  which  the  crossing  occurs. 

To  lay  the  foundation,  chapter  II  presents  a  review  of  the  Townes— Merritt 
effect,  quantum  beats  and  the  level  crossing  phenomenon.  The  original 
Townes— Merritt  analysis  used  a  semiclassical  approach,  which  is  not  compatible 
with  the  descriotion  of  level  crossings;  however,  the  results  are  useful  as  a 
comparison  to  a  quantum  electrodynamic  (QED)  model.  The  descriptions  of 
quantum  beats  and  level  crossings  provide  the  initial  mathematical  framework 
which  will  be  used  later  in  applying  the  Townes-Merritt  effect. 

Chapter  III  revisits  the  Townes-Merritt  effect  in  the  context  of  QED.  The 
Heisenberg  operator  formalism  is  used  to  describe  the  atomic  dynamics. 
Additionally,  the  calculations  use  the  adiabatic  approximation  to  describe  the  atom 
in  the  nonresonant  RF  field.  The  results  are  compared  to  a  dressed  states  model  to 
confirm  the  approach.  The  overall  purpose  is  to  understand  the  mechanisms 
involved  in  the  Townes-Merritt  sidebands  in  order  to  apply  them  to  the  cases  of 
quantum  beats  and  level  crossings. 

The  three  level  atom  is  pressed  in  chapter  IV.  Calculations  for  the 
instantaneous  energy  states  for  an  atom  in  an  RF  field  are  performed,  following  the 
adiabatic  approach,  for  both  linear  and  quadratic  Stark  shifts.  These  then  are  used 
in  descriptions  of  the  signals  resulting  from  quantum  beats  and  level  crossings. 
Expressions  describing  the  effect  of  the  RF  field  are  developed,  which  form  the  basis 
of  using  the  level  crossing  method  as  a  spectroscopic  technique. 

A  final  step,  in  chapter  V,  is  to  investigate  the  prospects  for  an  experimental 
test  of  this  method.  Suitable  atomic  transitions,  RF  field  strength  and  frequency, 
and  required  detection  capability  will  be  investigated.  The  intent  here  is  to  take  an 
initial  look  at  the  feasibility  of  conducting  this  type  of  experiment. 
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II.  Background 


Townes-Merritt  Effect 

As  a  preliminary,  a  review  of  t  he  original  Townes-Merritt  model  is  required 
since  it  will  provide  a  comparison  with  any  QED  approach.  Consider  an  absorption 
experiment  in  which  the  sample  is  placed  Li  a  rapidly  varying  nonresonant  field  and 
probed  with  a  second  resonant  beam.  The  original  calculations  used  semiclassical 
time-dependent  perturbation  theory  to  describe  transitions  from  this  dynamic  Stark 
effect  [2,7].  The  wave  equation  for  an  atom  in  the  nonresonant  field  is 

[Ho-zi-focos  J4,]  |^>  (2.1) 

where  H0  is  the  unperturbed  atomic  Hamiltonian;  p  is  the  dipole  operator;  £0  is  the 
nonresonant  field  amplitude;  and  v  is  the  field  frequency. 

Using  the  essential  states  approximation  with  two  unperturbed  states 
represented  by  1 1  )  and  |2  ),  the  solution  for  a  perturbed  state  can  be  written  as 

i  rl 

-If  f(t')dt' 

|0i  )  =  [a|l  >  +  b|2  >  Je  hJ0  (2.2) 

where,  assuming  the  perturbation  is  small,  a  a  1  and  b  is  small.  The  term  f(t)  is 
approximately  the  energy  of  the  unperturbed  state  Ei.  For  a  second  order  Stark 
effect,  the  wave  function  is  found  by  substituting  equation  (2.2)  into  equation  (2.1): 

j 0i  )  =  [a|  1  )  +  b]2  )]  expj—  ^(Ejt  +  +  ^-^-sin  2*4 )j  (2.3) 

where  AEi  is  the  static  Stark  shift, 
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In  the  absorption  experiment,  transitions  will  occur  between  |  ip j>  and 
another  st^ie  represented  by  |  $j).  If  the  interaction  is  H'(t)  =  //-ficos  uK, 
perturbation  theory  indicates  that  the  first  order  transition  probability  [8]  is 

t  2 

W’(t)|2  =  p  /  <  >dt'  (2.4) 

1  -00 

The  intensity  is  then  dependent  on  the  matrix  element: 

<V^|H'|^)  « /*3i£  exp|i(a»i3t  +  +  ^psin  2*4  ~  wt)J  (2.5) 

Significant  transitions  will  occur  only  when  the  exponential  term  is  "slowly 
varying,"  that  is  when  its  argument  is  approximately  zero.  The  exponential  can  be 

simplified  using  the  relation  e**28111  £  Jn^e**11^  [9]  yielding 

n  =  -<» 

<*|  H'  |  *>  «  w,£exp{.'(Wl,+  [^]e’2n‘/l  (2.6) 

Equation  (2.6)  describes  transitions  which  will  occur  for 

w  —  <^i3  +  — "b  ^nt*,  n  =  0,±1,±2,  .  .  .  (2-7) 

with  relative  intensity  Jn  .  This  is  the  case  for  an  atom  which  exhibits  a 
quadratic  Stark  effect.  The  calculation  describes  absorption  sidebands  resulting 
from  the  presence  of  a  nonresonant  electric  field.  The  position  of  these  sidebands  is 
primarily  dependent  on  the  field  frequency.  These  predictions  agree  with 
experimental  results  [7]  within  experimental  error.  For  the  case  of  a  linear  Stark 
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effect,  transitions  occur  for  u;  =  0^3  ±  ni/,  with  relative  intensity  Jl  >  where 
40-3  =  [2], 

Discovery  of  this  effect  has  led  to  promising  applications.  Arimondo  and 
Glorieux  [10]  and  Rackley  and  Butcher  [11]  performed  absorption  experiments  using 
a  nonresonant  RF  field  to  tune  an  atomic  transition  into  coincidence  with  a  laser. 

In  effect  they  adjusted  the  RF  frequency  until  the  fixed  laser  frequency  matched  one 
of  the  sidebands.  In  a  second  type  of  application,  Skatrud  and  DeLucia  [12] 
reported  on  using  an  RF  field  to  Stark  tune  a  laser.  Here  the  sidebands  of  an  atomic 
transition  were  used  for  the  lasing  transition  and,  hence,  the  laser  was  tuned  by 
changing  the  frequency  of  the  RF  field. 

In  each  of  these  examples,  application  of  the  Townes— Merritt  effect 
demonstrated  some  significant  advantages  over  a  static  or  DC  Stark  shift.  Because 
of  its  frequency  dependence,  as  per  equation  (2.7),  the  Townes— Merritt  effect  is  not 
limited  in  the  shift  magnitude  due  to  breakdown  voltage,  thus  allowing  a  wider 
range.  Since  the  DC  Stark  shift  depends  on  the  electric  field  amplitude,  a  uniform 
field  is  important  to  reduce  broadening.  The  Townes— Merritt  effect  does  not  rely 
on  field  uniformity,  as  the  intensity  only  determines  relative  intensities  of  the 
sidebands.  Finally,  the  measurement  or  control  of  field  frequency  is  substantially 
more  accurate  than  that  of  field  intensity.  This  allows  greater  experimental 
accuracy.  These  advantages  lead  to  the  proposed  application  of  the  Townes- 
Merritt  effect,  the  use  of  level  crossings  to  determine  fine/hyperfine  structure. 

Level  Crossings  and  Quantum  Beats 

Level  crossings  and  quantum  beats  are  related  quantum  phenomenon  that 
arise  from  interference  between  two  atomic  transitions  which  have  the  same  final 
energy  state  [13]  as  shown  in  Figure  4.  Level  crossings  occur  when  the  two  initial 


8 


Figure  4.  Three— Level  Atomic  System.  Transitions  can  occur  from 

either  excited  state  to  the  ground  state 

states  1 1  )  and  1 2  )  coincide,  resulting  in  a  change  in  the  polarization  and  angular 
dependence  of  the  fluorescence  pattern.  Quantum  beats  can  occur  when  the  initial 
states  are  slightly  different  in  energy;  the  result  is  a  modulation  of  the  total  radiated 
power  at  a  frequency  equivalent  to  the  energy  difference.  In  both  cases,  the  initial 
states  must  be  excited  coherently  to  achieve  the  interference  effect. 

In  level  crossing  experiments,  one  observes  the  resonant  fluorescence  of  the 
atom  as  a  function  of  an  external  field.  A  change  in  the  intensity  or  polarization  at 
a  specific  observation  angle  signifies  the  occurrence  of  a  level  crossing,  see  Figure  5 
[14].  The  crossing,  or  degeneracy,  is  due  to  Stark  or  Zeeman  shifts  in  the  original 
atomic  energy  states  resulting  from  an  external  static  field.  [6,15]  The  Hanle  effect 
is  a  special  case  of  a  zero  field  level  crossing.  The  primary  application  of  the  level 
crossing  technique  is  the  determination  of  original  energy  separation  based  on  the 
field  strength  at  which  crossing  occurs.  The  sensitivity  of  this  method  is 
comparable  to  the  natural  line  width  of  the  excited  states.  Several  experiments 
have  reported  using  this  method  [16,17,18]  to  measure  various  parameters, 
establishing  it  as  a  high  resolution  spectroscopic  technique.  One  of  the  drawbacks  is 
that  the  static  field  must  be  relatively  strong  and  uniform. 
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Figure  5.  Level  Crossing  Experiment  [6l.  Fluorescence  intensity  as  a 

function  of  magnetic  field  ana  the  atomic  energy  level  diagram. 

The  quantum  beat  phenomenon  occurs  when  there  is  a  slight  difference  in  the 
energy  of  the  initial  excited  states.  When  excited  coherently,  the  transitions  from 
these  states  interfere  resulting  in  a  modulation  of  the  radiated  power  at  the  "beat" 
frequency.  [191  Figure  6  shows  a  representative  example  of  a  beat  signal  in  a 
fluorescence  experiment.  [20]  A  principal  use  of  the  quantum  beat  effect  is  to 
investigate  very  slight  splittings  in  atomic  levels  [21,22]  where  other  spectroscopic 
techniques  do  not  have  sufficient  resolution. 

To  understand  the  theory  describing  level  crossings  and  quantum  beats, 
consider  the  simple  three  level  model  shown  in  Figure  4  in  which  two  excited  states 
can  decay  to  a  common  ground  state.  The  hypothetical  experiment  involves 
coherently  exciting  the  two  upper  states  and  observing  the  resonant  fluorescence. 
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CHAMMCl  NUW tn 

Figure  6.  Quantum  Beats  in  Zeeman  states  of  Cd:  5  3Pi,  [20].  Beats 
result  in  a  modulation  of  the  exponential  decay  curve. 

fluorescence.  The  Poynting  vector  describes  the  radiation  from  the  atom;  in  QED 
form  this  is 

S  =  ^{E<->x  B<*>  E<+>)  (2.8) 

where  the  superscript  plus  and  minus  refer  to  the  positive  and  negative  frequency 
parts  of  the  E  and  B  fields.  [23]  An  ideal  detector  measures  the  expectation  value  of 
this  expression  which  reduces  to 

S  =  <  E<-> -E<+>  >  (2.9) 

This  is  consistent  with  the  photo  ionization  process  which  takes  place  in  the 
detector,  described  by  Glauber.  [24] 
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From  dipole  radiation  theory  [25],  the  E  field  can  be  expressed  as 


= 


"(±) 
T  x  UK  ' 

c2r3 


(2.10) 


where  p  is  the  dipole  acceleration  operatoi.  From  Cook  [3],  an  expression  for 
radiation  intensity,  in  terms  of  the  dipole  operator,  results  from  substitution  of 
equation  (2.10)  into  equation  (2.9): 


'-‘(-I  ”<♦)  "(-)  V<+) 

s  =  r  <  //i  -pi,_  /ij_  j/ij.  )  (2.11) 

using  the  summation  notation  for  the  vector  components,  equation  (2.11)  gives  a 

means  of  predicting  the  fluorescence  radiation  at  any  point  by  the  expectation 

V  (  - )  V  (  ♦ ) 

value  of  the  second  time  derivative  of  the  dipole  operators,  ( fij  •/. i)  ).  This 
forms  the  basis  for  both  the  level  crossings  and  quantum  beat  descriptions.  In  level 
crossings,  the  item  of  interest  is  the  radiation  intensity  in  a  particular  direction  as  a 
function  of  the  energy  difference  in  the  excited  states.  For  quantum  beats,  the  total 
power  radiated  as  a  function  of  time  is  of  concern. 

First  consider  the  case  of  quantum  beats.  The  total  power  radiated  is  equal 
to  the  integral  of  equation  (2.11)  over  a  sphere  of  radius  r: 


P  = 


■<-)  -<+) 

4(  ft  ) 

3c=» 


(2.12) 


Now  to  find  an  expression  for  the  dipole  moment,  let  the  energy  eigenstates  of  the 
atom  be  defined  by 
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H|n)  =  E„|n> 


(2.13) 


These  states  are  orthogonal  and  complete.  Then  a  convenient  operator  basis  set  is 
(Tnm  =  |n)(m|,  such  that  in  the  Heisenberg  picture, 

kt)  =  (2H) 

D  ,  m 

where  Ham  =  (n|/i|m)  are  the  dipole  transition  moments.  This  can  be  separated 
into  the  separate  negative  and  positive  frequency  components  as 


-  <  -) 

/* 


n>m 


(2.15a) 


-<♦) 


Xu  <r  eI&/nn|t  =  \  u*  a*  eta;pq 

**nniffnme  pqc 

m>n  p>q 


(2.15b) 


Taking  the  second  derivative  with  respect  to  time  and  multiplying, 


”(-)  -<+) 

li  n 


\  \  V2  <JLu  -U  <r  e^^iua-^pq)1 
2^  “'nnr'pqPnm  ^pq^nm^pq^ 

n>m  p>q 


(2.16) 


The  last  term  in  this  equation  shows  the  presence  of  the  beat  frequency,  u>nm— Wpq- 
At  this  point  it  is  worth  noting  that  beats  will  not  occur  for  all  combinations 
of  transitions,  as  semiclassical  theory  would  predict.  Only  those  pairs  which  have  a 
common  final  state  exhibit  the  beats.  To  see  this,  consider  the  a  terms, 

ln><ml  [|p)<ql]f  =  !n)(mlq)(p|  =  ^np  (217) 
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which  shows  that  the  final  states,  |m  )  and  |q  )  must  be  the  same. 

equation  (2.16)  is  then  simplified  by  eliminating  the  sum  over  q: 


-<-)  -(♦) 

ft  • ft 


(2.18) 


n,  p>m 


where  u>„p  =  u>nra  —  u>pm  is  the  beat  frequency.  In  taking  the  expectation  value  of 
this  sum  in  the  state  j  ip  )  =  E  am  |  n  )  ,  we  find 


<  O'np  >  =  <^|n><p  |iii>  =  apa„  =  />po 


(2.19) 


which  is  the  density  matrix  in  the  Heisenberg  picture.  Hence,  the  expression  for  the 
total  power  becomes 


X  wnm  v 

n,  p>m 


2  U  U  pH  e*“'npt 

pm  J*pm  ^pn  w 


(2.20) 


As  an  example,  consider  a  three-level  atomic  system  as  shown  in  Figure  4. 
With  the  ground  state  defined  by  m=0,  there  are  tour  combinations  of  n,  p.  The 
radiated  power  is  then  given  by 


P  =  Molgiol  _pt  +  li4oi^2oLp2 


4u%(\u{ol  *  H  W 12t  .  *  H  iu> 21 1 

- —  PlO’^20  P2l€  ^  +  P20‘/^10  Pl2 e 


(2.21) 


This  clearly  demonstrates  the  presence  of  a  beat  frequency  U21  ,  provided  the  upper 
states  are  coherently  excited,  p!>i  ^  0. 


For  level  crossings,  specific  characteristics  of  the  radiation  field  as  a  function 
of  the  energy  level  separation  are  important.  Generally,  a  change  in  the  fluorescent 
intensity  at  a  given  observation  point  indicates  the  level  crossing.  To  develop  an 
expression  for  this,  rewrite  equation  (2.11)  as 


c  (r26i i-xjxi)  /  “< *>  * (*\ 

S  =  1  <  >“  -W  > 


(2.22) 


Using  the  expression  for  the  dipole  product  in  equation  (2.18) 


S  =  ^~2rc4^  X  wpm  Km- 4* 

n,  p>m 


(2.23) 


Once  again  consider  the  three  level  atomic  system  shown  in  figure  4.  As 
before  the  summation  yields  four  terms.  Assuming  £Uio  ~  ^20  *  wo  , 


#*10  #*io  P1  +  4)  tio  P2 

+  Mio4o  /»2l(t)  +  /4oMio  pM 


(2.24) 


This  describes  the  fluorescence  as  a  function  of  position.  It  depends  on  the  dipole 
moments  and  the  density  matrix. 

To  gain  further  insight,  consider  a  relatively  simple  case  of  Zeeman  splitting 
in  the  2p  states  of  hydrogen,  specifically  the  m  =  ±1  states.  The  dipole  moments 
for  these  states  can  be  written  [27] 


Mw  =  u0(x  +  *y) 

ho  =  M0(x  ~  iy)  (2.25) 
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If  the  spontaneous  emission  rate  is  A  for  both  states,  the  density  matrix 
elements  associated  with  the  upper  levels  can  be  written  as: 

P.  =  P2  =  5^Q2e_Al  (2.26a) 

f>i2  =  *>2i  =  -2^Q2e'At~At  (2.2Gb) 

where  Q  is  the  time  integral  of  the  exciting  pulse  envelope  divided  by  h  and  A  is  the 
frequency  corresponding  to  the  energy  difference  between  the  excited  states  [28]. 
Using  these  equations  in  equation  (2.24), 

S  =  [1  +  cos2*  +  sin2*  cos(2<^>-At)]  (2.27) 

where  *  and  <(>  define  the  observation  angle,  in  spherical  coordinates. 

A  detector  placed  at  a  given  angle  will  see  an  energy  pulse  whose  total 
energy  is  defined  by  integrating  equation  (2.27)  over  time: 


J  4  jrc^r2 


-1  _+_  cos J  +  A^2^(Acos  2<f>  +  Asin  2<f>) 


(2.28) 


For  the  direction  *  =  ^  ,  <f>  =  Q  this  becomes 


n  _  Wo/ToQ2 


1  I  A 

A  +  A2  +  A2 


(2.21)) 


Figure  7  shows  equation  (2.29)  as  a  function  of  the  energy  difference.  There  is  a 
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Figure  7.  Fluorescent  Intensity  at  Level  Crossing.  Crossing  occurs  at  A  =  0. 

significant  change  in  the  observed  fluorescence  at  A  =  0.  In  addition,  this  has  a 
narrow  width,  on  the  order  of  the  natural  linewidth  (FWHM=2A),  making  this 
technique  potentially  useful  for  high  resolution  spectroscopy. 

The  mathematical  descriptions  of  quantum  beats  and  level  crossings  in  this 
section  provide  the  necessary  framework  to  address  the  effects  of  a  nonresonant  RF 
field  on  the  atomic  fluorescence.  Equation  (2.20)  for  quantum  beats  and  equation 
(2.28)  for  level  crossings  describe  signals  as  a  function  of  the  energy  difference 
between  the  excited  states.  Similar  expressions  will  be  developed  for  the  atom  in  an 
RF  field,  except  that  the  energy  difference  will  be  a  function  of  the  field  frequency 
due  to  the  Townes-Merritt  effect.  The  next  section  examines  the  Townes-Merritt 
effect  using  the  Heisenberg  operator  formalism  to  describe  the  atomic  dynamics. 
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III.  QED  Description  of  the  Townes— Merritt  Effect 

To  understand  the  impact  of  the  Townes— Merritt  effect  on  level  crossings 
and  quantum  beats  requires  a  reformulation  of  the  Townes— Merritt  description  in 
terms  of  quantum  electrodynamics  (QED).  It  is  convenient  to  use  the  Heisenberg 
operator  formalism  in  developing  a  model  to  describe  the  effect.  Initially,  a  two 
level  atom  in  a  resonant  field  is  examined  to  gain  insight  in  using  the  formalism. 
Then  the  case  of  the  atom  in  a  nonresonant  field  is  discussed.  Here  the  adiabatic 
approximation  is  introduced,  which  will  be  central  in  the  description  of  level 
crossings  and  quantum  beats.  Combining  the  nonresonant  field  with  a  near 
resonant  field,  using  the  adiabatic  approximation,  results  in  a  description  of  the 
Townes— Merritt  effect.  Two  cases  will  be  considered:  a  linear  Stark  shift  of  a  single 
excited  state  and  a  quadratic  Stark  shift.  A  QED  model  is  then  developed,  and  the 
results  compared  to  the  adiabatic  approach.  This  validates  the  adiabatic  approach 
as  a  means  to  describe  the  atomic  dynamics  in  a  nonresonant  field.  As  a  final  step, 
the  dressed  states  model  is  used  to  describe  the  Townes— Merritt  effect.  The  results 
further  corroborate  the  adiabatic  solution  which  will  be  used  in  level  crossings  and 
quantum  beats. 

Heisenberg  Operator  Formalism. 

In  the  Heisenberg  operator  formalism  the  set  of  operators  defined  by 

<^ij  =  I »  > ( j I  (:}1) 

forms  a  complete  operator  basis,  where  11  j  i  )  =  E{  |  i  ).  An  arbitrary  operator  A 
can  then  be  described  by 
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A  =  (3.2) 

*.  j 

where  Ajj  h  (  i|  A(0)|  j  ).  In  the  Heisenberg  picture,  then,  the  o  operators  obey  the 
equation  of  motion 


aij  ~  tfr  [aij’Hl  (3-3) 

a  H 

with  the  initial  condition  a  j  j(0)  =  |  i  )  (  j  j .  Thus,  the  time  evolution  of  any 

A  U 

operator  can  be  determined  by  solving  for  o’jj(t),  and  the  expectation  value  of  the 
operator,  which  represents  a  physical  observable,  can  be  found  by  taking  the 

A  U 

expectation  values  of  o  Vy 

The  a  operators  exhibit  useful  mathematical  properties.  They  are  a  closed 
set  under  multiplication 


<Tij<Tnm 


(3.4) 


which  is  useful  in  evaluating  equations,  such  as  commutation  relations.  The 
commutation  rule  for  these  operators  is 


(3.5) 


In  addition,  the  Hamiltonian  for  an  unperturbed  atom  can  be  represented  by 


II 


o 


(3.6) 
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Figure  8.  Two  Level  Atom 

*  ii 

The  time  evolution  of  ox:  is  then  given  by 


E2  —  iho>0 


Ei=  -tfK 


*ij(0  = 


(3.7) 


where  w-  =  . 

Consider  now  the  two  level  atom  with  energy  eigenstates  shown  in  figure  8. 
Define  the  a  operator  basis  in  the  following  way: 

*  =  U  X  2| 

^3=  1 2  ) <  2 1  —  |1  >(  1|  (3.8) 

These  two  operators  contain  all  the  information  for  the  two  level  atom.  The 
unperturbed  atomic  Hamiltonian  can  be  represented  as  H0  =  £hu>0a3 .  In  the 
electric  dipole  approximation,  the  interaction  Hamiltonian  for  the  atom  in  the 
electric  field  is  given  by 


H,  =  -  ji-E(t) 


(3.9) 
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The  electric  dipole  approximation  is  valid  when  the  E  field  spatial  variation  is 
negligible  over  the  extent  of  the  atom  [8],  which  is  the  case  here  as  the  RF  and 
optical  wavelengths  are  large  compared  to  the  atomic  diameter.  The  form  of  /i  and 
E  will  define  specific  cases  for  the  interaction.  Initially  consider  the  atomic  levels  to 
have  definite  parity  and  treat  the  E  field  classically.  The  dipole  moment  is  then 


t*  = 


0 


]c  =  n{a 


+ 


at)f 


(3.10) 


where  e  is  the  unit  orientation  vector.  Let  the  electric  field  be  classically  prescribed 
in  the  form  E  =  efcos  u t  ,  which  is  oriented  in  the  same  direction  as  the  atomic 
dipole.  For  a  strong  field  this  is  reasonable  as  will  be  shown  in  the  QED  treatment. 
The  total  Hamiltonian  is  then 


*  i  «.+ 

H  =  ^hu)0<T3  ~  i£(o  +  a  )cos  cnt  (3.11) 

Using  the  Heisenberg  equation  of  motion  and  the  commutation  relations  for  the  a 
operators, 


a  —  —iuca—  iQ. cos  ut  <r3  (3.12a) 

1  »  -  + 

a3  =  — i2flcos  uA  (a  —  a  )  (3.12b) 


where  H  ^  is  the  Rabi  frequency. 

If  the  field  frequency  is  close  to  the  transition  frequency,  a  u0,  a  near 
resonance  condition  exists.  The  first  term  on  the  right  hand  side  of  equation  (3.12a) 
is  the  dominant  term,  so  an  approximate  trial  solution  can  be  written 
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C  —  Ik'r.t, 

a  =  Se  0 


(3.13) 


where  S  is  slowly  varying  in  time.  Substitution  into  equations  (3.12)  yields 

S  =  -^(e*^  e“ttft)e,w‘>t  a3  (3.14a) 

a3  =  e~iu^)  (5e-*^ot-  5te*wot)  (3.14b) 

Since  u>  a  u0,  the  terms  with  frequency  of  ~2 u0  can  be  ignored,  as  they  are  rapidly 
oscillating.  These  do  not  contribute  significantly  to  the  solution  since  they  average 
to  zero  over  a  short  time  period,  much  less  than  the  period  of  an  observation.  This 
is  commonly  referred  to  as  the  rotating  wave  approximation.  Thus, 

s  =  4fe-,At  u, 

ir3  =  -if!(ie'At-5te_‘At) 

where  A  =  u  -  uj0  is  the  detuning. 

Of  particular  interest  is  the  on  resonance  case,  when  A  =  0.  Equations  (3.15) 
can  then  be  reduced  to 


(3.15a) 

(3.15b) 


<73  =  —  fi203  (3.16) 

A  solution  for  the  inversion  operator,  a3  ,  can  easily  be  found,  after  taking  the 
expectation  value  of  equation  (3.16): 

<r3(t)  =  <73(0)  sin(flt  +  <j>)  (317) 
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This  shows  the  Rabi  oscillations  at  frequency  ft.  The  Rabi  frequency  depends 
directly  on  the  field  intensity,  and  a  higher  frequency  indicates  greater  absorption. 

In  order  to  more  precisely  define  the  absorption  rate,  spontaneous  emission 
needs  to  be  considered.  From  a  phenomenological  argument  [29],  spontaneous 
emission  can  be  included  as  additional  terms  in  equations  (3.12) 

*  *  a* 

a  =  —  iu0a  —  tftcos  o^t  <r3  —  (3.18a) 

<r3  =  — i2ftcos  wt  (<r  —  a  )  —  A(<t3  +  1)  (3.18b) 

Justification  for  these  terms  will  be  shown  later  in  the  QED  section  where  the 
electric  field  is  quantized. 

Following  an  approach  similar  to  the  previous  case,  let  a  =  Se  lu^  .  The 
equations  of  motion  then  become 

£  =  iA£-ifttf3-y<7  (3.19a) 

<r3  =  — ift(5  —  5f)  —  A(a3  +  1)  (3.19b) 


For  steady  state  conditions,  S  =  a  =  0,  the  solution  for  the  expectation  value  of  a3 
is 


A2  +  4A2 

—  —  A2  +  T7F  +  STP 


(3.20) 


The  probability  of  being  in  the  upper  state  is  then 

p2  =  S*1  +  =  hJ  + 


(3.21) 
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Figure  9.  Absorption  Cross-Section.  The  lineshape  is  Lorentzian,  with  a 
weak  field  width  (FWHM)  of  A. 

This  is  a  Lorentzian  function  of  width  (FWHM)  {A2+2ft2.  The  term  2 H2 
represents  power  broadening.  For  weak  field  conditions,  A2  >  202,  the  probability 
of  being  in  the  upper  state  is  directly  proportional  to  H2.  The  absorption 
cross-section  a A  is  equal  to  the  rate  of  spontaneous  emission,  AP2,  divided  by  the 
incident  photon  flux  density  .  I  is  the  light  intensity  which  can  be  expressed  as  I 
Thus, 

-  _  hwAP2  _  8  twAm2  _ _ _ 

a  T"^  -  ficfA2  +  4A2  +  2fi2J  (3-22) 

where  H2  =  .  Figure  9  shows  the  absorption  cross-section  as  a  function  of  the 

detuning  A.  The  A  coefficient  can  be  expressed  as  A  =  ^  ,  so  for  a  weak  field 
the  peak  is  o~  =  6?r — n  . 

0  Wo2 

Adiabatic  Approximation 

The  Townes-Merritt  effect  occurs  when  there  are  two  electric  fields  present, 
one  nonresonant  which  will  be  referred  to  as  the  RF  field  and  one  resonant  probe 
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field  with  a  frequency  near  the  transition  frequency,  u  »  u0.  The  RF  field  frequency 
v  is  much  less  than  the  transition  frequency.  The  Hamiltonian  is  then 


H  =  K0-iiEv-iiK  (3.23) 

To  find  a  solution  describing  the  atomic  dynamics,  first  consider  only  the  effect  of 
the  RF  field  on  the  atom.  Since  it  is  far  off  resonance  and  is  slowly  varying  as 
compared  to  the  transition  frequency,  the  adiabatic  approximation  is  useful  [8,30]. 
Here  the  instantaneous  energy  eigenstates  form  a  convenient  basis  set: 

HA(t)|n(t)>=E„(t)|n(t)>  (3.24) 

where  H  (t)  =  H  —  ft- E„.  The  physical  interpretation  of  this  approach  is  that  as 
the  RF  field  changes,  the  eigenstates  evolve  adiabatically  with  the  field.  The  state 
|n(0)  )  at  time  zero,  with  field  strength  zero,  will  evolve  to  a  state  |n(t)  )  at  time  t 
that  is  Stark  shifted  by  the  field  E(t).  Since  the  RF  frequency  is  far  off  resonance,  it 
is  reasonable  to  assume  no  transitions  occur.  Additionally,  these  eigenstates  remain 
orthonormal  throughout.  The  total  wave  function  is  written  as 


Itft)  >  =  ^ajn(t)  ) 

n 


(3.25) 


This  approach  is  valid  so  long  as  an,  which  represents  transitions,  is  small. 

As  stated,  the  adiabatic  approximation  is  useful  for  cases  in  which  the 
Hamiltonian  changes  slowly  in  time,  in  comparison  to  the  atomic  dynamics.  This  is 
just  the  case  for  the  atom  in  an  RF  field  when  the  frequency  v  is  much  less  than  the 
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transition  frequency  o/jj.  In  the  adiabatic  treatment,  the  energy  eigenstates  and 
eigenvalues  are  smooth  functions  of  time  which  change  in  a  continuous  manner  with 
the  Hamiltonian  satisfying  the  'time  independent'  wave  equation  (3.24). 

Of  particular  interest  is  the  Heisenberg  equation  of  motion  for  a5j-  =  |i(t)  )( 
j(t)  | .  In  the  adiabatic  approach,  the  equation  of  motion  will  be  approximated  as 


(3.26) 


The  additional  terms  will  be  be  small  and  rapidly  oscillating,  such  that  they  do  not 
contribute  to  the  evolution  of  the  states  from  which  physical  transitions  occur. 

Conversion  to  the  Heisenberg  picture  is  defined  by  a  unitary  transformation 
U(t)  which  satisfies 


l*t)>  =  U(t)llK0)>  (3.27) 

where  |  tp  )  is  the  solution  to  the  Schrodinger  wave  equation.  Substituting  equation 
(3.27)  into  the  wave  equation 

ft^W)>  =  H(t)U(t)|«0)>  (3.28) 

From  this,  the  time  dependence  of  U  can  be  seen: 

6(t)  =  H(t)U(t)  (3.29) 


An  operator  is  defined  in  the  Heisenberg  picture  by 
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Ah  =  Ut  As  U 


(3.30) 


Now  Take  the  derivative  and  use  equation  (3.29)  to  get 


a"  =  Uf  AS  U  +  Uf  As  U  +  U*  As  U 


IriH  ,  3A’ 

=  riA  ’Hl+5r 


(3.31) 


nAH  *  +  A  g  *  c 

where  the  last  term  is  defined  as  =  U  A  U  ,  which  arises  if  A  has  an  explicit 
time  dependence.  This  then  is  the  equation  of  motion  for  any  operator  in  the 
Heisenberg  picture. 

*  ii 

For  o  j j,  using  equation  (3.24),  the  equation  of  motion  becomes 


<r”j  =  tWij(t)<T?j  + 


(3.32) 


where  oi^t)  =  .  An  expression  for  the  last  term  can  be  developed  by 

first  finding 

^=|i><j|  +  |i><j| 


dfrii  _  ii  \/  ;i  m  i;  \/  ;i  (3.33) 


Using  the  identity  operator 


^»i=£|k><k]i><j|  +  £|i><j|m><m| 

k  m 

=  X<kl‘>  ^kj  +  X<jlm>  (3‘34) 

k  m 
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Now  an  expression  for  the  inner  products  is  found  by  differentiating  equation  (3.24) 


Then  multiply  by  (  k| ,  where  k# 


'  '  '  fu^i  k 


For  k  =  i,  (  i|i  )  =  0  since  (  i|i  )  =  1.  Thus, 


,  +  Y  <  iffi"  >  ft 

**  riwj 


m#j  ^ j  m 


If  the  RF  field  is  written  as  E(t)  =  £cos  i/t,  then 


dH  •  c7E(t)  *  r  • 

=  —  /*•  —Qi  =  P'S  I'sm  iA 


W 


The  equation  of  motion  then  becomes 


*?i  =  S'W^i  +Xi%irfwin  ^  *!j  +Xa^rein 

k#i  ,k  a*j  Jm 


(3.35) 


(3.36) 


(3.37) 


(3.38) 


(3.39) 


The  first  term  on  the  right  hand  side  is  the  dominant  term,  which  is  seen  by 
comparing  magnitudes  of  the  coefficients. 


(Hi*' 


^ik 


<  Kl 


(3.40) 
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where  Qki  =  is  the  Rabi  frequency,  and  is  typically  less  than  v.  This  is  the 

condition  of  validity  for  the  adiabatic  approximation  [30].  As  a  result,  an 

»  1/ 

approximate  solution  for  ^(t)  can  be  written 

^j(t)  =  SH{ t)  e^i1  (3.41) 

where  is  a  slowly  varying  function  and  uP-  is  the  unperturbed  transition 
frequency.  Substituting  this  into  equation  (3.39) 

K  =  0ifsin  14  V“k,t 

k  i  i 

+  (3.42) 

mi]  Jm 

Notice  that  all  terms  in  the  summations  have  frequencies  of  order  uki  ±  v.  Since 

v  <  u)ki,  these  terms  are  all  rapidly  oscillating  and  can  be  neglected.  The  result  is 

Q  h 

that  the  terms  arising  from  — ^  can  be  ignored  in  the  equation  of  motion,  which 

A  U  «  |j 

can  then  be  approximated  as  al}  s  iu;ij(t)  al}  for  the  atom  in  an  RF  field.  The 
solution  takes  the  form 

-h  -  H  »/iij(t')  dt' 

<4j(t)  =  »",(0)  e  J°  (3.43) 

This  result  will  be  most  useful  in  describing  the  atomic  dynamics  in  a  nonresonant 
field.  When  a  near  resonant  field  is  added,  the  Townes-Merritt  effect  results.  For 
level  crossings  and  quantum  beats  spontaneous  emission  is  considered.  But  in  all 
cases,  evaluation  of  equation  (3.43)  for  a  specific  set  of  energy  states  is  central. 
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The  adiabatic  approximation  defined  a  set  of  energy  eigenstates  for  the  atom 
in  an  RF  field.  Application  of  a  near  resonant  field  will  now  cause  transitions  to 
occur  between  the  energy  states.  These  transitions  will  occur  not  only  at  the 
fundamental  frequency,  but  also  at  frequencies  shifted  by  integer  multiples  of  the 
RF  frequency;  this  is  the  Townes— Merritt  effect.  The  analysis  here  is  very  similar 
to  that  of  an  atom  in  a  resonant  field,  except  in  this  case  the  energy  states  are  the 
instantaneous  eigenstates  from  the  adiabatic  approach. 

Consider  the  two  level  system  depicted  in  figure  8.  An  exact  solution  exists 
for  the  energy  states  when  an  RF  field  is  applied  [27]: 

cos  0\  \  sin  0 

l»  >  =  [_(„,]  |2>=[cos«j  <M4> 

where  tan  20  =  —  ^cos  vi  .  This  assumes  a  quadratic  Stark  shift,  the  linear  case 
U)0 

will  be  covered  later.  Note  that  since  usually  <  u0  ,  values  of  0  will  typically  be 
small,  the  new  states  will  be  approximately  the  same  as  the  original.  The 
instantaneous  energy  levels  are 

E±  =  ±  rj1  [J*0  +  4172cosVt]*  (3.45) 


where  the  plus  and  minus  refer  to  the  upper  and  lower  states  respectively.  The 
instantaneous  frequency  difference  is  then  given  by 


K  +  402cos2t't]’ 


(3.46) 


In  the  Heisenberg  operator  formalism  the  a  operators  will  be  defined  by  these 


adiabatic  states  such  that, 


*ij(0  E  l«(t)  X  j(t)|  (3.47) 

and  their  time  evolution  will  in  turn  define  the  atomic  dynamics. 

The  Hamiltonian  of  equation  (3.23)  can  be  rewritten  as  H  =  HA  +  , 

where  the  interaction  with  the  resonant  field  is  given  by  =  - /i£  cos  uA.  The 
Heisenberg  equation  of  motion  for  the  a  operators  is  then 


'u  =  i  M*  +  H  J  =  m  J  (3-48) 

where  the  h  superscripts  for  the  Heisenberg  picture  have  been  dropped.  Now 
consider  the  case  of  the  two  level  atom  with  a  operator  defined  by  equations  (3.8): 

a  =  -iV(t)  a  -  tfi^cos  wt  a3  (3.49a) 

<r3  =  —i2U^ cos  uA  (a  —  a)  (3.49b) 

These  equations  are  very  similar  to  the  case  of  only  the  resonant  interaction  except 
for  the  time  dependence  of  uf.  The  solution  follows  much  the  same  approach. 
Assume 

ifl'(V)  dt1 

<r(t)  =  5(t)  e  (3.50) 

This  is  reasonable  since,  in  equation  (3.49a),  —iui'(i)a  is  the  major  contributing 
term.  Substituting  the  assumed  solution  into  the  equation  of  motion  and  discarding 


31 


the  rapidly  oscillating  terms  (rotating  wave  approximation),  one  obtains 


*  o  -«/  A(t')  dt*  . 

s=-q»eJ  0  a, 

f-  *  r  A(f)  dt’  -  +  —if  A(t')  dt' 
*3  =  0  0 


(3.51a) 

(3.51b) 


where  A(t)  e  u-  [w2  4-  4fl2cos2i4,]’.  First  consider  the  exponential  terms, 


i f  A(t')  dt'  if  u-  +  4n2cos2J4']’dtl 

•'a  —  p  "  O 


e  0  =  e 


(3.52) 


Since  flv  «  u0  , 


[uP0  +  4ft2C0S2K]’  »  a/0  +  C0S2t4,  + 


(3.53) 


Ignoring  higher  order  terms  leads  to 


X 


r  {Uo  +  202  cos’*')dt'  =  w0t  +  jgt  +  sin  2* 


(3.54) 


An  exponential  term  with  a  sine  argument  can  be  expanded  as  a  Bessel  function 
series: 


02 

-t-'—VTfiin  2 lA  ^2, 


—  It; — Sill  n2 

e  ^  -  IJn(^) 


Ql  s  -:2n/4, 
e 


(3.55) 


n  =  -oo 


The  exponential  then  becomes 
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(3.56) 


tJ^A(t')  dt' 


i(u>  —  u0 
e 


—  2w)t 


Now  consider  an  on  resonance  case  where 


Then, 


Similarly, 


n? 


uKu)m  =  u0  +  -f  +  2nw 
0 


XA(t')dt'= 

-./o  A(t')  c-e(m-„)K 


The  equations  of  motion  become 


(3.57) 


(3.58) 


(3.59) 


-  o  *5"  o 


(3.60a) 

(3.60b) 


For  the  case  of  interest  the  RF  frequency  is  much  greater  than  the  Rabi 
frequency,  v  >  ,  so  that  power  broadening  does  not  cause  individual  sidebands  to 

overlap.  For  a  similar  reason,  the  RF  field  frequency  must  be  greater  than  the 
natural  linewidth,  v  »  A.  Since  the  operators  evolve  at  the  frequency  ,  the  only 
contributing  factors  in  the  expansion  series  is  for  m  =  n  ;  the  others  rapidly  average 
to  zero.  This  is  the  equivalent  of  the  rotating  wave  approximation.  Thus,  the 
equations  of  motion  become 
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(3.61a) 


=  (3.61b) 

These  are  similar  to  the  previous  case  without  the  RF  field,  equations  (3.16),  except 
for  the  Bessel  function  coefficients.  With  a  zero  RF  field  strength,  they  revert  to 
the  earlier  form.  The  solution  for  the  inversion  operator  <r3  is 

=  ff3(0)  sinlJ.ljSL)  +  <j>]  (3.62) 

Again,  this  shows  Rabi  oscillations  at  frequency  reflecting  the  relative 

absorption  strength  of  the  sidebands.  The  result  agrees  with  the  original 
Townes-Merritt  calculation  for  a  quadratic  Stark  shift.  The  absorption  occurs  at 

Q2  2  f)2 

frequencies  um  =  w0  +  — "  4-  2mv  ,  with  relative  intensity  • 

To  include  the  effects  of  spontaneous  emission,  consider  a  more  general  case 
of  the  equations  of  motion  where  the  incident  radiation  is  slightly  off  resonance: 


(3.63a) 
(3.63b) 

where  again  A(t)  a  w  —  [u%  +  4f)2cos2/4,]’.  Now  let  w  =  w0  +  ^  +  2mv  +  Aw  , 

where  a ui  is  small,  describing  the  condition  where  w  is  slightly  detuned  from  the  m 
sideband.  The  exponential  terms  can  be  expanded  similar  to  the  Following  the 


5  =  -fif*  e  Jo  -  A 


S 


<t3  —  —in 


r-  *  r  A(t’)  dt’  -+  -if  A(t')  dt’ 
—S'e  Jo 


-  A(a3+1) 
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development  of  equation  (3.59),  which  results  in 


±* rtA(f)  df 

e 


±tAii>t 

=  e 


e±i2(m— n)i4, 


(3.64) 


Elimination  of  all  terms  in  the  expansion  except  for  m  =  n,  results  in 


s  =  oz  ^ -§S  (3.65a) 

*3  =  -iSVU^)  t**4*  -sV^]  -A(ff3+1)  (3.65b) 

*  *  I  *.  jL 

To  remove  the  remaining  exponential  terms,  let  S  =  S  e  .  The  equations  of 
motion  for  the  expectation  values  are  then 


S  -  i&uS  -  J,n(5^p)  <^3  -  ^  ^  (3.66a) 

*3  =  — aWsfj)  [«'- $'*]  -A(ff,+1)  (3.66b) 

The  steady  state  solution  for  the  inversion  operator  is 

_  _  A2  +  4AU>2  to  ro\ 

a3  -  A2  T  TO  +  2^2 J2  (3-67) 


This  result  is  similar  to  that  without  the  RF  field,  equation  (3.20).  The  only 

2  Q2 

difference  is  the  Bessel  function  term  Jm(<r-f/7)  which  serves  to  determine  the 

Li  U/qI' 

relative  strengths  of  the  sidebands.  For  fl„  =  0  ,  Jm(0)  ~  ^Om  and  the  equation 
reverts  to  the  non— RF  case. 
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Figure  10. 


Absorption  cross-section  for  an  atom  in  a  nonresonant 
field  of  frequency  v. 


The  absorption  spectrum  consists  of  a  series  of  peaks  separated  by  2v.  Near 
each  peak  the  probability  of  being  in  the  upper  state  is  given  by 


P  Vi  i  T  1 

P2  -  s(l+<73)  -  A2  +  4a^  + 


(3.68) 


where  Acj  is  the  local  detuning.  From  equation  (3.22),  the  absorption  cross-section 
is 


8  iru)Av2J  m 

°a  “  hc(A2  +  4Syr  +  2Jjyaj 


(3.69) 


n2 

Each  peak  is  a  Lorentzian  function  centered  at  u  =  u0  +  v  4-  2mv ,  having  a 


linewidth  (FWHM)  \  A2+2fl2J£  .  For  a  weak  field,  Q2  <  A2,  the  individual 
absorption  peaks  are  <r0J2  ,  where  o0  is  the  peak  absorption  in  the  non  RF  case. 
The  spectrum  is  shown  in  figure  10,  which  shows  the  expected  absorption  sidebands. 
This  agrees  with  the  original  Townes-Merrittt  analysis  [2,7], 
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The  case  for  the  linear  Stark  shift  is  very  similar  to  the  quadratic  case, 
except  the  calculation  for  the  instantaneous  energy  states  is  much  simpler,  as  will  be 
shown.  Recall,  the  Hamiltonian  for  the  atom  in  an  RF  field  is 


HA  =  -  M*  K 


(3.70) 


For  a  linear  shift,  the  interaction  term  can  be  written 


P-EV  = 


/hi  0 

P  1*22] 


£cos  vi 


(3.71) 


To  determine  the  energy  levels,  let  H  |n  )  =  En|n  )  ,  which  results  in 


Ej  =  -  ?jha;0  -  f/iucos  ui 
E2  =  -  <f/i22cos  lA 


(3.72a) 

(3.72b) 


The  instantaneous  transition  frequency  is  then 


a;'(t)  =  u0  —  (lv cos  lA 


(3.73) 


where  . 

As  in  the  previous  case,  when  a  probe  beam  is  turned  on,  the  equations  of 
motion  for  the  Heisenberg  operators  become 

a  =  — iV(t)  a  —  tf^cos  A  (73  (3.74a) 

<r3  =  -i2Qjcos  A  (a  -  c?)  (3.74b) 
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The  solution  is  similar  to  the  previous  case.  Define  A(t)  =  to  -  u>0  4-  Ovcos  i/t  , 
then 


if\(  t')dt' 
e 


*  ,o  ,  i(w-w(>  +  ni/)t 

1  w) e 

n=  -oo 


(3.75) 


Let  uj  ~  u>0  —  mi/ ,  and  consider  only  terms  where  n  =  m  : 

(3.76a) 
(3.76b) 


s  J„(S»)  *3 


This  reduces  to  o3  =  — D^J^(^)  a3  ,  and  the  solution  can  be  expressed  as 


a, 


*3  =  *j(°)  sin(Jm(l7,/)  ^ 


(3.77) 


This  differs  from  the  result  in  the  quadratic  case  by  the  Bessel  function 
argument  and  that  the  transitions  occur  for  every  integer  multiple  of  the  RF 
frequency  u  with  no  static  shift,  instead  of  just  the  even  harmonics.  The  expression 
for  steady  state  absorption  in  the  linear  case  is  the  same  as  equation  (3.69)  except 
again  the  Bessel  function  term  has  a  different  argument.  The  relative  strength  of 
the  absorption  sidebands  is  given  by  J^(^)  • 

The  descriptions  of  the  Townes-Merritt  effect  derived  using  the  Heisenberg 
operator  formalism  and  the  adiabatic  approximation  agree  exactly  with  previous 
published  approaches  and  experimental  results  [2,26].  This  gives  a  great  deal  of 
credibility  to  the  approach,  which  will  be  used  to  describe  level  crossings  and 
quantum  beats.  But  before  moving  on  to  those  subject,  it  is  important  to  consider  a 
quantum  electrodynamic  analysis  of  the  Townes-Merritt  effect. 
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luantum  Electrodvnamic  Treatment 

In  the  previous  Section,  spontaneous  emission  was  included  based  on  a 
phenomenological  argument.  A  more  rigorous  approach  involves  quantization  of  the 
radiation  field.  This  QED  approach  gives  stronger  justification  for  the  spontaneous 
emission  terms,  and  provides  a  more  accurate  representation  of  the  atom-field 
system.  Recall,  that  level  crossings  and  quantum  beats  result  from  QED  effects.  It 
will  be  shown,  however,  that  for  a  strong,  classical-like  field  the  QED  description  of 
the  Townes— Merritt  effect  agrees  with  the  Heisenberg  operator  approach  presented. 

The  Hamiltonian  for  the  atom-^ield  system  in  QED  is 

H  =  jH(t3  +  ^Tfu^a^a^  -  ii-E  (3.78) 

The  first  term  is  the  Hamiltonian  for  the  isolated  atom.  The  second  represents  the 

~  j- 

electromagnetic  field  where  akX  and  akX  are  the  creation  and  destruction  operators. 
The  third  term  is  the  interaction  Hamiltonian  in  which  the  electric  field  is  now 
quantized.  In  the  Townes— Merritt  effect,  the  electric  field  consists  of  the 
nonresonant  RF  field  and  a  near  resonant  probe  beam: 

E  =  E„  +  Ew  (3.79) 

First  consider  the  atom  and  the  RF  field.  Following  the  adiabatic  approach, 
define  a  set  of  basis  states  which  satisfy  the  equation 

HA|n(l))=E„(t)|n(t)>  (3.80) 

where  HA  =  jhai0<r3  —  ft- E„  ,  as  before.  These  states  are  the  instantaneous  energy 
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eigenstates  of  the  atom  in  the  RF  field,  which  adiabatically  follow  the  relatively 
slowly  changing  field.  In  the  Townes-Merritt  effect,  the  RF  field  is  sufficiently 
strong  that  ran  be  treated.  Thus,  expressions  for  the  energy  states  and  values 
would  agree  with  those  previously  developed.  Also  recall,  in  the  adiabatic 
approximation,  the  equation  of  motion  for  aVl  is 

°\]  =  ft  Kj>HAl  =  *u>ij(t)<Tij  (3.81) 

Here  the  h  superscript  for  the  Heisenberg  picture  has  been  dropped  for  simplicity. 
Adding  the  probe  field,  ,  the  equation  of  motion  becomes 

*ij  =  K’H]  =  ^  EJ  (3.82) 

-  - 

The  operator  aV]  commutes  with  the  a  and  a  operators  so  the  electromagnetic  field 

term  of  the  Hamiltonian  has  no  impact  on  the  equation  of  motion.  To  derive  an 

.  .  „  + 

expansion  for  the  remaining  commutator  relation,  recall  p  =  p(a  +  o  )  .  Also  for 
the  two  level  system,  use 


a  =  |1  )<  2| 

*3  =  |2  ><  2|  -  (1  ><  1| 


as  the  a  basis  set.  The  equations  of  motion  are  then 


a  =  -t'u/(t)  a-^/i-E^  <r3 
=  °)  VK 


(3.83a) 

(3.83b) 


(3.84a) 

(3.84b) 
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where  u/(t)  =  .  The  total  field  commutes  with  the  atomic  operators 

since  it  is  a  function  of  a  and  a*. 

The  electric  field  can  be  separated  into  free  field  anu  source  field  terms: 

E*  =  E0+Efi  -  (3-85) 

The  free  field  term  represents  incident  electromagnetic  radiation.  The  source  field 
results  from  the  radiation  reaction  field  of  the  atom.  These  terms  can  be  further 
separated  into  positive  and  negative  frequency  components,  after  Glauber  [24]: 


E=  E(+)  +  E(') 


(3.86) 


This  is  easily  seen  in  the  quantized  expression  for  the  free  field: 


k,  X 


(3.87) 


where  the  first  term  with  a  is  E^+)  and  the  second  term  E),"' .  Before  discussing  the 
source  field,  there  are  some  relations  between  the  a  and  a^  operators  and  field 
coherent  states  to  cover,  which  are  important  in  evaluating  the  equations  of  motion. 

A  coherent  state  provides  the  closest  quantum  mechanical  representation  of  a 
classical  electromagnetic  wave.  It  is  a  good  description  of  laser  radiation,  which  is 
the  case  of  interest  here.  The  coherent  state  is  a  linear  combination  of  number 
states  given  by 


o )  =  e” 


a] 


n4tfnT 


n  ) 


(3.88) 
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where  a  is  any  complex  number.  Unlike  for  the  number  states  |n),  the  expectation 
value  of  the  E  and  B  fields  are  non-zero  for  coherent  states.  Coherent  states  are 


also  described  as  minimum  uncertainty  states,  for  the  product  of  the  position  and 
momentum  uncertainties  is  the  minimum  ^  allowed  by  the  Heisenberg  uncertainty 
principle.  This  type  of  state  is  then  the  closest  description  of  a  strong  classical  field, 
l  he  coherent  sittie  ib  <ui  eigen  veCi.C-  of  the  uestructicn  operator  a: 

a|  a  )  =  a  )  (3.89) 

The  field  coherent  state,  |  a  )p  ,  is  defined  as  the  product  of  individual  coherent 
states: 

l«>F=  K>K>  •  •  •  (3-90) 

where  the  indices  represent  all  possible  k  and  A.  Relations  for  field  coherent  states 
are 

akxl«)F  =  <*kxl«>F  (3-91a) 

F(  ttlakX  =  akXF(  ttl  (3.91b) 

With  these  expressions,  it  can  be  shown  that 

F<  =*o_,F<  a| 

where 

k,  X  ^ 

£x 


(3.92a) 

(3.92b) 

(3.93a) 

(3.93b) 


Equation  (3.92)  will  be  useful  in  developing  equations  of  motion  for  the  a  operators. 

The  source  field,  Eg ,  occurs  due  to  the  presence  of  the  atom.  It  is  also 
referred  to  as  the  self  field  or  radiation  reaction  field.  It  arises  from  the  need  to 
compensate  for  energy  lost  from  the  system  due  to  radiation.  Classically,  as  the 
atom  radiates  energy,  the  energy  of  the  atomic  dipole  must  decrease  accordingly. 
Quantum  mechanically,  the  expectation  value  of  the  electric  field  due  to  the  atomic 
dipole  decreases.  From  classical  electrodynamics  and  the  correspondence  principle, 
a  QED  expression  for  the  source  field  is  [31] 


Es  -  _  3^3  V  (a  ~  °)  (3.94) 

This  agrees  with  the  QED  derivation  [29].  In  the  equation  of  motion  for  a,  equation 
(3.84a),  the  primary  term  on  the  right  hand  side  is  -u;'(t)<r  since  -J  >  .  Thus, 

the  third  derivative  can  be  approximated  as 

<7  ss  tu'3(t)  a  (3.95) 

The  source  field  expression  then  becomes 

E(  =  E«*>  +  E>-’  =  (3.96) 

which  will  be  used  in  developing  the  equations  of  motion. 

The  equations  of  motion  are  put  in  the  normal  form,  where  the  a^  operators 
are  to  the  left  and  the  a  operators  are  to  the  right,  by  separating  the  electric  field 
into  positive  and  negative  frequency  portions: 
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a 

*3 


=  -iu/'(t)  <r- ^  ft-  E^erj  -f  E(  +) | 


=  -2* 


(/i-  E(  "*  (<7  —  O^)  +  (<7  —  E(  +) 


(3.97a) 

(3.97b) 


This  is  possible  since  the  positive  and  negative  frequency  portions  individually 
commute  with  the  atomic  operators.  Each  field  term  can  then  be  split  into  fiee  field 
and  source  field  terms.  Using  equation  (3.96)  and  the  a  operator  product  rule 


a  =  —  *u7'(t)  a  —  ^  /i-E(  _)  <r3  +  E(  +)  1  —  ^{a  -  er^) 


^3  ~ 


-t 


*tx 


p-E —  o  )+(a—<r  )/i-E<+)  —  A(<t3+  I) 


(3.98a) 

(3.98b) 
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where  A  =  is  the  Einstein  A  coefficient,  representing  the  rate  of 

spontaneous  emission.  This  leaves  only  the  free  field  terms. 

Since  the  atom-field  system  is  fully  quantized,  the  state  of  the  system  can  be 
represented  by  the  product  of  the  atomic  state  and  the  field  state: 


l*>=  l^A>k>F  (3-99) 

where  |  a  )p  is  a  field  coherent  state.  Taking  the  expectation  values  and  using  the 
eigenvalue  relations  of  equations  (3.92),  the  equations  of  motion  are 

<7  =  ~iuf(t)  cr-£/i-£0(t)a3-j(a-  a  )  (3.1™a) 

ff3  -  -2jr/i-50(t)[a-  a*}  -  A(a3  +  1)  (3.100b) 

where  £0(t)  =  £(0*>  +  £(n~'  ,  represents  the  incident  proi/i  beam  and  can  be  written 
as  £0(t)  =  £0cos  uA.  In  this  case,  the  equations  of  motion  are  equivalent  *o  those  in 


44 


the  previous  section  using  the  adiabatic  approach.  To  demonstrate  this,  let 
a  =  <S(t)  exp  j— tj'^Xt')  dt1  J  ,  where  <£(t)  is  a  slowly  varying  function.  Since  u  a  (J  . 

the  rotating  wave  approximation  is  valid,  and  the  equations  become 


S  -- 


0  —  i  f  A(t')  dt 

e 


A  e 

*3- j* 


A(<r3+1) 


(3.101a) 

(3.101b) 


where  and  A(t)  =  u  —  w'(t)  .  These  are  equivalent  to  equations  (3.63) 

in  the  adiabatic  approach;  hence,  the  results  are  the  same. 

For  a  quadratic  Stark  shift,  the  instantaneous  transition  frequency  is 
u/(t)  =  [u%  4-  fi2cos2i/t]’.  Transitions  then  occur  a  frequencies  -  u0  +  ^  + 

U/Q 

02 

2mv  ,  where  m  is  an  integer,  with  relative  intensity  J2(^— 1 v~)  ■  The  instantaneous 

4  uJqi/ 

transition  frequency  for  a  linear  shift  is  u>'(t)  =  u0  -  flvcos  iA  ,  and  transitions 
occur  at  u;m  =  u0  +  mi/  with  relative  intensity  J2^)  . 

The  results  of  this  more  rigorous  description  of  the  Townes-Merritt  effect 
agree  exactly  with  the  previous  method  and  with  the  original  Townes-Merritt 
experiment.  This  adds  greatly  to  the  confidence  of  the  model  using  the  Heisenberg 
operator  formalism  and  the  adiabatic  approximation  to  describe  the  effect  of  the 
nonresonant,  RF  field.  This  adiabatic  approach  will  be  used  to  model  a  three  level 
atom  in  level  crossings.  But  first,  a  third  method  using  dressed  states  to  describe 
the  Townes-Merritt  effect  will  be  examined  to  provide  further  insights  to  the 
physical  processes. 
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Figure  1 1  Dressed  States  for  a  Two  Level  Atom.  Each  atomic  state  is 
combined  with  a  photon  ladder  representing  the  possible 
number  of  photons  in  the  field.  The  combined  energy  states 
are  represented  by  dashed  lines. 

Dressed  States  Model 

Dressed  states  are  useful  in  visualizing  the  transitions  which  can  occur  in  the 
Townes— Merritt  effect,  although  the  mathematics  is  cumbersome.  The  purpose  here 
is  to  describe,  with  the  dressed  states  model,  transitions  which  give  rise  to  the 
absorption  sidebands,  and  thus,  further  validate  the  adiabatic  approach.  Cases  for 
both  the  linear  and  quadratic  Stark  shift  will  be  examined. 

A  given  dressed  state  is  a  combination  of  the  atomic  state  and  the  photon 
number  state.  The  energy  level  diagram  consists  of  the  atomic  energy  states 
superimposed  with  a  ladder  of  the  photon  states.  Figure  11  shows  a  two  level  atom 
in  a  nonresonant  field.  The  state  |  l,m  )  represents  the  ground  atomic  state  with  m 
photons  present,  and  is  equal  to  the  product  of  the  two  states,  1 1  )a|m  )  .  When 
a  probe  beam  is  applied,  transitions  may  occur  between  the  ground  state  ladder  and 
the  excited  state  ladder,  say  between  |  l,m  )  and  |2,n  )  .  Allowed  transitions 
depend  on  the  actual  energy  states.  If  only  the  original  states  are  considered,  the 
dipole  matrix  would  be 
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(3.202) 


/*im2n  =  (  l,m|^|2,n  >  =  /ij2<  m|n  )  =  p, 2Smn 


This  would  prescribe  that  only  transitions  where  n  =  m  could  occur.  However,  the 
actual  energy  states  are  linear  combinations  of  the  unperturbed  states, 

ImO  =  XCiinm^’m  >  (3‘103) 

j .  m 


so  transitions  where  n  ^  m  may  occur.  These  states  are  indie*.'  ed  by  dashed  lines  in 
figure  11,  and  are  slightly  shifted  in  energy.  The  transition  strength  is  proportional 
to  the  square  of  the  dipole  moment  and  is  then  dependent  on  the  coeffi  lent  Cjjnni  . 
The  transition  frequency  is  a  function  of  the  new  energy  levels.  The  new  states 
satisfy 

H|MT>  =  Ein|MT)  (3.104) 


and  thus,  are  dependent  on  the  Hamiltonian.  Two  cases  will  be  considered  here,  the 
linear  Stark  shift  of  the  excited  state  and  a  quadratic  Stark  shift. 

For  the  linear  case,  the  Hamiltonian  can  be  written  as 


H  =  hu>0</cr  4-  hv& a  -  /*0-E|2)(2j(a  +  a^)  (3.105) 


where  |E|  = 


hi/ 


TTa V 


=  £v/2[m  and  m  is  the  average  number  of  photons  [29].  The 


ground  state  is  unchanged,  so  |  l,m  )  =  j  l,m  )  .  The  new  excited  state  ladder  is 
a  linear  combination  of  the  |2,m  )  states  only,  given  by 


l2’n  >  =  XC"J2’m  > 


(3.106) 
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Note  that  the  coefficient  Cnm  should  tend  towards  £nm  as  the  field  strength  goes  to 
zero.  Substituting  this  into  equation  (3.104) 


]£Cnm(M>+  to™)|2,m)  -  ^Cnm  M0‘E[fm|2,m-l>  +  ^m+l|2,m+l>] 

m  m 

=  E2n£cnm|2,m)  .  (3.107) 

m 

Multiplying  by  (  2,m'  |  results  in 

C„.(K+  Tu/m)  “  Cnm-i  +  Cnm+1  ft0‘Efm  =  E2nCnm  (3.108) 

f 

Let  E2n  =  hu>0  +  nh^ ,  and  assume  \  m+1  «  (m  w  |m  which  is  valid  for  high 
photon  occupation  states.  Equation  (3.108)  then  becomes 

Cn„a2Trk=cn„.1  +  cn„.1  (3.109) 

where  f!  =  =  s\  as  before.  This  recursion  relation  is  equivalent  to  one 

for  Bessel  functions  (9],  where  CnD)  =  .  Thus  the  energy  eigenstates  are 

|5nr>  =  l2.m  >  <311°) 

with  energy 

E2n  =  hu0  +  nfw  (3.111) 

Note  as  £  tends  toward  zero,  Jm-n(^)  becomes  6nm  . 

Transitions  from  the  atomic  ground  state  to  the  excited  state  are  a  function 
of  the  dipole  matrix  element.  For  an  atom  initially  in  the  j  1  m  )  state,  the 
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transition  probability  to  the  |2,n  >  state  is  proportional  to 


=  K  l,mjp|2,n  >|2 

=l^o!2|XJt  "(§,')  <m|k) 

k 

-  U2  j2  Aa 
ro  Jm*nV  p  ) 


(3.112) 


The  transition  would  occur  at  frequency  u>  =  u0  +  (n— m)fu/ ,  with  relative 

2  ft 

intensity  •  This  agrees  with  the  result  in  the  adiabatic  approach;  here  the 

difference  term,  n— m,  represents  the  n— m  sideband.  For  the  total  transition 
strength,  the  individual  moments  are  summed  over  the  distribution  in  the  ground 
state  ladder,  but  the  relative  strengths  of  the  sidebands  would  remain  the  same. 

For  the  quadratic  case,  the  situation  is  a  little  more  complex  since  the  energy 
eigenstates  are  functions  of  both  original  atomic  states  summed  over  the  photon 
states: 

\h*~)=  £cjnkmlkim>  (3-113) 

k,  di 


The  Hamiltonian  for  a  many  level  atom  in  the  RF  field  is 

H=  +  ^ta-^ann/inn*E[a  +  at]  (3.114) 

n  n  ,  m 

Substituting  equations  (3.113)  and  (3.114)  into  (3.104) 

XCjnkm(K+  mfu/)|k,m)  -  ]^Cjnkmppk-E[Jm|p,m-l>  +  jm+T  |p,m+l>] 

k >  m  p,  k,  tn 

=  Ejn^Cjnkmlk>m>  (3I15) 

k ,  m 
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Multiplying  by  (  k',m'  j  results  in  the  following  recursion  relation: 

(K  +  mfu/—  Ejn]Cjnkm  =  ^  /*kP  *  E[ jm+TCj  n  pm+ ,  +  jHCjnpnrl]  (3.116) 

p 

This  is  still  a  general  relation;  the  linear  Stark  shift  of  the  previous  example  occurs 
when  ^  is  non-zero  only  for  k  =  p  =  2.  For  the  quadratic  shift,  pkp  =  0  for  k  =  p. 

Consider  initially  the  upper  state  of  the  two  level  atom,  j  =  2  .  Let 
Ej,,  =  hu>0  +  nhv  +  en  .  The  index  k  can  take  on  values  of  1  or  2.  For  k  =  2, 

[(m-n)fu/—  en]C2n2m  =  /ij,  •  E[^m+TC2nln)+1  +  JmC^.J  (3.117) 

For  k  =  1, 

[-h^0+(m-n)fu/ -  fn]C2nlm  =  ft2-E[fm+lC2n2m+1  +  (iHC2n2m.1]  (3.118) 

Then  substitute  equation  (3.118)  into  equation  (3.117): 


r  /•„  ,  lr-  _  l..  17*  l  2  r(rn-+-l  )C2n  2  m  +  (  m+2  )(m+l  )  C2r 

[(m— n)7u^—  cn]C2n2m  —  l/*2i  EJ  w  +  er 


mC2n2m  +  <1  m(m-l  )C2n2ni-2 
hw0  -  (m-n-l)lu/  -  en 


(3.119) 


Again  assume  a  highly  populated  photon  state,  m  a  m+1  a  m  ,  and  let  tn  = 

,  where  11  5  .  Noting  that  hu0  >tw, cn  ,  the  recursion  relation 

simplifies  to 

^  u°  C2n2tn  =  C2n2m*2  +  C2n2m.2  (3.120) 

This  relates  alternate  coefficients  in  the  photon  ladder  of  the  excited  state.  Since 
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the  perturbation  is  relatively  small,  the  largest  coefficient  is  for  m  =  n.  In  fact, 
when  the  RF  field  goes  to  zero  this  is  the  only  coefficient  in  the  series.  In  the 
presence  of  an  RF  field,  the  other  non-zero  coefficients  could  only  be  the  even 
series,  or  where  m  =  n  ±  2k;  k  is  an  integer.  This  result  agrees  with  the  experiment; 
absorption  sidebands  exist  only  for  the  even  harmonics  of  v. 

Equation  (3.120)  is  a  Bessel  function  recursion  relation  where  C2n2m  = 

Q2 

JJs-n(^jt,~) ,  imposing  the  condition  m  =  n  ±  2k  ,  for  the  even  series.  The  energy 
states  can  then  be  written 


i2^~> = i2’n+2k  >  <3-m> 

k  0 

with  energy 

E;„=K  +  nfu/+^  (3.122) 

o 

This  ignores  terms  involving  the  ground  atomic  state  since  the  coefficients  C2nIm 
are  smaller  by  a  factor  of  ,  as  can  be  seen  in  equation  (3.118).  The  last  term  in 
the  expression  for  the  energy  is  the  static  Stark  shift,  which  agrees  with  previous 
theory.  Similar  expressions  for  the  ground  state  ladder  can  be  developed: 

ith> = XJi(^ 1 1,m~2j }  (3123) 

j 

=  (3.124) 

Now  consider  a  transition  from  the  j  l,m  )  state  to  the  |2,n  )  state  as  a 
result  of  the  probe  beam.  The  strength  of  this  transition  is  directly  proportional  to 
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The  last  line  results  from  the  relation  m  —  2j  =  n  +  2k  .  Notice  also  that  for  this 
to  be  satisfied,  the  indices  m  and  n  must  differ  by  a  multiple  of  2,  m  -  n  =  2p  . 

This  becomes  a  selection  rule  for  the  transition.  The  sum  of  the  Bessel  function 
products  can  be  further  reduced  by  the  relation  E  Jk(z)Jv_k(z)  =  Jv(2z)  [9];  thus, 

l/hm2nl  ~  $2  (3.126) 

!•  O 

_  ry 

Transitions  can  then  occur  at  frequencies  <j\-  %  +  +  2pt> ,  where  p  =  is 

2  Q2 

an  integer.  The  relative  strengths  are  given  by  . 

The  result  here  agrees  with  previous  approaches  in  describing  the 
Townes— Merritt  effect,  which  provides  an  added  degree  of  confidence  in  the 
adiabatic  approximation  for  defining  the  energy  states  of  an  atom  in  a  nonresonant 
field.  The  dressed  states  model  provides  a  convenient  picture  for  understanding  the 
physical  transitions;  however,  its  complexity  warrants  use  of  the  adiabatic  approach 
in  describing  the  three  level  atom  for  level  crossings  and  quantum  beats.  The 
dressed  states  picture  raises  one  other  issue,  the  question  of  quantum  beats  in  a  two 
level  atom  due  to  interfetence  between  transitions  from  adjacent  atom— field  states, 
which  will  now  be  discussed. 
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(uantum  Beats  in  a  Two  Level  Atom 

Quantum  beats  have  not  been  observed  in  a  two  level  atom,  nor  are  they 
predicted  in  the  adiabatic  model.  In  the  dressed  states  picture,  the  possibility  arises 
if  adjacent  photon  states  in  the  excited  atomic  state  ladder  transition  to  a  common 
state  in  the  ground  state  ladder.  These  type  of  transitions,  n,p  photons  to  m 
photons,  are  allowed  due  to  the  mixing  of  states,  as  shown  in  the  previous  section. 
This  is,  however,  a  simplistic  view  of  the  transitions  from  excited  state  to  ground 
state  since  it  ignores  the  total  contribution  from  all  states  in  the  ladder.  It  will  be 
shown  that  when  all  possible  transitions  are  considered,  there  is  no  beating  effect, 
although  the  sidebands  are  still  present. 

The  quantum  beat  phenomenon  occurs  when  two  excited  states  can 
transition  to  a  common  ground  state.  Interference  between  the  two  possible 
trant.  .ions  can  produce  a  beat  signal.  In  the  dressed  states  picture,  it  is  possible  in 
a  twc  level  atom  fcr  two  different  excited  states  to  transition  to  the  same  ground 
state,  if  one  considers  the  various  atom-field  states.  To  demonstrate  this,  consider 
the  simplest  example,  where  only  the  upper  state  exhibits  a  linear  Stark  shift. 

Recall  |l,n  )  =  |  l,n  )  ,  and  from  equation  (3.110), 

1^0  =  Jjk.n(§-)  |2,k  >  (3.127) 

k 

The  mixing  of  states  allows  this  situation  to  exist.  Figure  12  depicts  relaxation 
from  various  levels  in  the  excited  state  ladder  to  a  single  level  in  the  ground  state 
ladder,  which  could  conceivably  produce  the  beat  signal  at  multiples  of  the  RF 
frequency.  To  fully  evaluate  the  potential  for  interference,  these  type  of  transitions 
must  be  summed  over  the  entire  excited  ladder. 
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Figure  12  Quantum  Beats  in  a  Two  Level  Atom.  Transitions  occur 
between  combined  atom-RF  field  states. 

The  quantum  beat  signal  is  given  by  equation  (2.20),  which  is  repeated  here, 
assuming  w10 


p  =  <3128) 

n  Tp?m 

The  dipole  matrix  elements  are  between  the  excited  state  (2^n  j  and  the  lower  state 
|  l,m).  If  the  field  is  considered  to  be  a  coherent  state,  then  it  can  be  assumed  that 
the  excited  states  of  interest  have  equal  initial  populations.  This  is  reasonable  since 
in  a  strong  field,  the  coefficients  of  the  coherent  state  change  little  between  adjacent 
states,  and  the  coefficients  of  the  energy  states  |2,n  )  drop  off  rapidly  as  k  differs 
from  n. 

The  dipole  matrix  can  be  written  as 

Mnm=  <2,k|pjl,m  > 

k 

=  /vWp,/)  (3129) 
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Upon  substitution  of  this  into  equation  (3.128),  the  power  radiated  is 


p  =  X  P)VJ„ 

n 


The  Bessel  function  series  reduces  to  E  Jm-nJp-n  =  £np  ;thus  , 


(3.130) 


(3.131) 


Since  P  is  constant,  there  are  no  quantum  beats  present  despite  the  existence 
of  transitions  on  the  various  harmonic  sidebands.  An  interpretation  of  this  is  the 
individual  beat  signals,  which  would  be  present  if  the  atom  were  initially  excited  to 
only  two  states  in  the  ladder,  become  washed  out  when  combined  with  all  other 
possible  transition  pairs.  Quantum  beats  can  still  occur  in  a  three  level  atom  as  will 
be  shown  later. 
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IV. 


Quantum  Beats  and  Level  Crossings 


Three  Level  Atom 

Level  crossings  and  quantum  beats  result  from  interference  between  the 
spontaneous  decay  of  two  coherently  excited  states  which  have  the  same  final  or 
ground  state.  To  treat  this  phenomenon,  the  simplest  atomic  model  is  that  of  a 
three  level  atom  with  the  two  excited  states  close  in  energy;  refer  to  figure  4.  A 
level  crossing  occurs  when  the  upper  states  coincide,  resulting  in  a  shift  of  the 
fluorescence  pattern  depending  on  the  original  transition  characteristics.  When  the 
energy  states  are  slightly  separated,  quantum  beats  may  result  which  modulate  the 
radiated  power  during  spontaneous  emission. 

When  the  atom  is  placed  in  a  nonresonant  RF  field,  additional  transitions 
can  occur,  as  discussed  in  the  section  on  the  Townes— Merritt  effect.  These 
transitions  give  rise  to  other  possible  level  crossings  or  quantum  beats,  which  is 
precisely  the  subject  of  this  study.  The  adiabatic  and  the  dressed  states  models, 
used  to  describe  the  Townes— Merriti  effect,  provide  the  mathematical  basis  and  an 
intuitive  understanding  of  the  physical  processes  occurring.  The  dressed  states 
model  gives  a  clear  picture  in  terms  of  identifying  transitions  which  fit  the  three 
level  model  for  quantum  beats  or  level  crossings.  Figure  13  shows  the  photon 
ladders  associated  with  each  atomic  state.  Here  a  level  croosing  condition  exists 
when  the  |  l,n  )  and  |2,p  )  states  coincide  and  can  each  decay  to  the  |0,m.)  state. 
This  of  course  represents  but  one  of  the  possible  combinations  in  the  dressed  states 
model.  To  arrive  at  the  total  solution,  all  possible  combinations  need  be  included. 
This  added  complexity  is  the  prime  motivation  for  using  the  adiabatic  approach  to 
mathematically  describe  the  level  crossing  and  quantum  beat  effects.  In  the 
description  of  the  Townes-Merritt  effect,  the  adiabatic  solution  agreed  with  the 
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Figure  13  Dressed  states  picture  of  a  three  level  atom  in  an  RF  field. 

Energy  levels  represent  an  atomic  states  combined  with  a 
photon  state  with  n  photons  present. 

dressed  states  solution  as  well  as  the  original  experimental  results,  which  gives  the 
approach  significant  credibility. 

The  primary  difference  between  the  Townes-Merritt  experiment  and  either 
level  crossings  or  quantum  beats  is  that  in  the  latter  there  is  no  second  applied  field 
to  Uiggei  «t  transition.  Both  level  crossings  and  quantum  beats  result  from 
spontaneous  emission  of  the  atom  after  it  has  been  excited  by  a  short  pulse; 
however,  the  energy  eigenstates  which  result  from  the  RF  field  are  the  same.  Thus, 
the  calculations  will  consider  spontaneous  emission  transitions  from  the  excited 
energy  eigenstates  to  the  ground  state  of  the  three  level  atom. 

The  first  requirement  is  to  determine  the  instantaneous  energy  eigenstates 
for  the  atom  in  the  RF  field.  Two  cases  will  be  considered:  a  linear  Stark  shift  of  a 
single  excited  state  and  a  quadratic  Stark  shift  of  all  three  states.  The  first  case  is 
relatively  simple  and  the  atomic  operators,  can  be  determined  exactly.  The 
second  requires  an  approximate  solution,  using  perturbation  theory.  These  solutions 
will  then  be  used  to  calculate  the  quantum  beat  and  level  crossing  effects. 
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Linear  Stark  Shift.  For  the  linear  Stark  shift  case,  two  methods  will  be  used 
to  calculate  the  atomic  operators,  aVj.  The  first  solves  the  time  independent 
Schrodinger  equation  for  the  instantaneous  energy  eigenstates,  then  uses  them  to 
construct  the  a  operators.  The  second  uses  the  unperturbed  atomic  states  as  a  basis 
and  solves  the  Heisenberg  equation  of  motion  for  the  a  operators  in  the  Heisenberg 
picture.  Both  methods  result  in  the  same  expression,  but  the  latter  can  be  easily 
adjusted  to  include  spontaneous  emission. 

The  Hamiltonian  for  the  atom  in  the  RF  field  is  given  by 

n 

where 

/<= y 

i>  j 

For  the  case  of  a  single  level  linear  Stark  shift,  the  only  non-zero  term  of  nl}  is 
/i,!  =  fie.  Let  the  RF  field  be  given  by  E  =  e£cos  lA.  Then  the  Hamiltonian 
becomes 

H  =  ^hu;n<7nn  -  Wlv<rncos  lA  (4.3) 

n 

The  instantaneous  energy  eigenstates  are  solutions  of 

H(t)|n(t)>  =  E„(t)|n(t)>  (4.4) 

which  yields  |n(t)  )  =  |n)  and 


(4.1) 

(4.2) 
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Eo  =  0 

(4.5a) 

Ej  =  Uoj1  -  hftvcosiA 

(4.5b) 

E2  *”  fttdg 

(4.5c) 

The  a  operators  are  then  defined  by  <7^  =  |i(t)  >(  j(t)|.  In  the  Heisenberg 
picture, 

i|j  =  (4-6) 

The  solution  of  this  equation  results  in 


(4.7) 


which  defines  the  time  evolution  of  the  a  operators. 

The  second  method  starts  by  defining  |  i  >  (  j  |  ,  where  j  i  )  represents 

an  unperturbed  atomic  state.  The  Heisenberg  equation  of  motion  gives 


=  +  *1<4*u]^C0S  V t 

n 

=  +  ^ji»li  -  ^n<4j]S!„COS  *4 


(4.8) 


The  solution  of  equation  (4.8)  is  the  same  as  the  previous  method.  For  example, 
consider 

<7,2  ~  icJ,2<r , 2  —  iflvCO S  lA  <7,2  (4.9) 
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The  solution  is 


^(t)  =  £tJ2(0)  e 


(4.10) 


where  u>l2(t)  =  un  —  H„cos  iA. 

These  solutions,  however,  have  not  considered  spontaneous  emission.  This 
can  be  included  by  adding  a  decay  term  to  equation  (4.9) 


<t12  =  i<J12<r i2  —  tnvcos  iA  al2  —  A<t,2  (4.11) 

Here  it  is  assumed  for  simplicity  that  the  decay  coefficients  for  the  two  transitions 
are  equal  to  A.  The  result  agrees  with  the  previous  treatment  of  spontaneous 
emission  in  the  section  on  the  Townes— Merritt  effect,  which  lead  to  equation 
(3.98a).  The  factor  of  2  difference  arises  from  the  fact  that  both  the  1 1)  and  1 2) 
states  are  excited  states.  The  solution  then  reads 


(T^nU)  =  <7^2(0)  e 


-At  dt' 


(4.12) 


This  will  be  useful  in  calculating  the  effects  of  quantum  beats  and  level  crossings. 


Quadratic  Stark  Shift.  For  a  quadratic  Stark  shift,  the  methods  used  in  the 
previous  section  to  find  the  energy  states  do  not  yield  simple  solutions.  In  fact,  the 
second  method  involving  the  Heisenberg  equation  of  motion  would  require  solving 
six  coupled  differential  equations.  The  primary  reason  for  this  is  the  form  of  the 
dipole  matrix,  which  can  be  written 
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0  /*01  /*02 
f*  =  #*10  0  0 

#20  0  0 


(4.13) 


This  assumes  that  transitions  between  the  ground  state  and  each  of  the  excited 
states  are  allowed,  but  not  between  excited  states.  Additionally,  each  state  is  of 
definite  parity  so  =  0. 

The  time  independent  Schrodinger  equation  can  be  rewritten  in  matrix  form, 
using  equation  (4.13),  to  read 

0  fljCOS  ut  fl2cos 

h  ii^cos  lA  u>x  0  |n  )=En|n  >  (4.14) 

yi-icos  iA,  0 

where  •  Ae  exact  solution  requires  the  characteristic  equation,  in  this 

f 

case  a  cubic  equation  of  En  . 

A  more  tractable  approach  is  to  use  time  independent  pertuibation  theory  to 
find  the  instantaneous  energy  states.  This  is  valid  for  small  perturbations,  when 
Q  €  v  .  To  second  order,  the  perturbed  states  are  given  by  [8] 


where  H'  =  —  ;i-K  .  For  the  three  level  atom,  equations  (4.15)  and  (4.16)  yield 
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(4.17a) 


|(T>=  |0>-^kos«4|l>-^osI4|2) 


* 


=  l>>+^s«|0>+2;^2cos«|2> 

(4.17b) 

=  I2>+^*I°>  +  ^?4C0S1/‘I1> 

(4.17c) 

E0  =  -  cos2* 

0  aq  u>2J 

(4.18a) 

E,  =  tw,  +  Ml'™2* 

(4.18b) 

E,  =  +  M1’C0)!M 

(4.18c) 

From  these  expressions,  the  a  operators  may  be  written.  Note  that  in  forming  <rv.  , 
Cl  €  to  ,  so, 

*ij  =  U  X  j|  «  |i  X  j|  (4-19) 


3a ? 


Thus,  the  term  can  be  ignored,  and  the  Heisenberg  equation  of  motion  is,  as 


before 


au  -  iwyCt)  ajj 


(4.20) 


Spontaneous  emission  is  again  included  by  adding  an  exponential  decay  term  in  the 
solution.  In  the  case  of  a)2 ,  the  result  is 


where 


—At  1 f  )  dt’1 

a12(t)  =  a12(0)  e  Ate  Jo 


u>,2(t)  =  Lun  +  - — ^2cos2tA 

Uf\U> 2 


(4.21) 


(4.22) 
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Equation  (4.21)  is  the  same  result  as  for  the  linear  case,  equation  (4.12), 
except  for  the  definition  of  the  instantaneous  transition  frequency.  These  equations 
will  be  central  to  calculating  level  crossing  and  quantum  beat  effects  in  the  next 
section.  It  is  worth  noting  that  the  effects  for  the  quadratic  shift  are  expected  to  be 
much  less  due  to  the  smaller  instantaneous  shifts;  of  order  ^  less.  Also  if  Q,  =  fl2, 
the  effects  will  be  negligible. 


Quantum  Beats  Calculation 

Before  addressing  level  crossings,  the  effect  of  the  Townes— Merritt  sidebands 
on  quantum  beats  will  be  discussed.  Specifically,  this  section  will  cover  the  presence 
of  quantum  beats  in  which  the  magnitude  and  frequency  are  dependent  on  the  RF 
field.  This  is  a  simpler  calculation,  primarily  since  it  considers  the  total  radiation 
from  the  atom  as  a  function  of  time  instead  of  radiation  in  a  specific  direction. 

Both  the  single  linear  Stark  shifted  level  and  the  quadratic  Stark  shift,  covered  in 
the  last  section,  will  be  considered.  The  calculation  follows  the  adiabatic  approach, 
but  for  the  simpler  linear  case,  the  dressed  states  model  is  also  used  as  a 
comparison. 

Quantum  beats  can  occur  when  the  atom  is  initially  excited  by  a  pulse  such 

that  there  is  a  finite  probability  of  being  in  either  upper  state.  In  this  way,  the  off 

diagonal  elements  of  the  initial  density  matrix  are  non-zero.  The  off  diagonal  terms 

are  defined  as  the  product  of  the  amplitude  coefficients  of  the  wave  function  for 

* 

each  of  the  excited  states,  for  example  pn  =  ata2  .  As  the  atom  decays  back  to  the 
ground  state  from  the  'mixed'  upper  state  via  spontaneous  emission,  there  is 
potential  interference  between  the  two  distinct  transitions.  The  interference  results 
in  a  beating  effect,  or  an  amplitude  modulation  of  the  radiated  power,  which  has  a 
frequency  directly  related  to  the  energy  separation  of  the  two  excited  states.  The 


Townes-Merritt  effect  causes  a  splitting  of  these  excited  states  raising  the 
possibility  of  other  beat  signals  dependent  on  the  RF  field  frequency.  The  following 
calculations  will  show  this  dependence. 

Equation  (2.12)  described  the  power  emitted  from  the  atom  as  a  function  of 
the  expectation  value  of  the  inner  product  of  the  positive  and  negative  frequency 
portions  of  the  dipole  acceleration  operator: 


=  4 

3c3 


(4.23) 


In  the  adiabatic  approach,  the  dipole  operator  can  be  represented  as 

(4-24) 

n  ,  m 


where  <rnm(t)  is  defined  by  the  adiabatic  energy  states.  Following  the  development 
in  chapter  II,  an  expression  for  the  inner  product,  similar  to  equation  (2.18),  results: 


=  Y  u>2  u2 

nmwpm 

n  ?p?m 


^nm'^pm 


(4.25) 


The  expression  for  the  power  radiated  then  becomes 

P  =  3F  S  <  anP(t)  >  (4’26) 

n  rp">m 

Considering  the  possible  values  for  the  indices  in  the  case  of  the  three  level  atom, 
the  index  m  must  be  zero,  and  (n,p)  can  be  (1,1),  (1,2),  (2,1),  or  (2,2).  The  first 
and  last  possibilities  are  trivial  since 
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<  i„(t)  >  =  P,(t)  =  P,(0)  e  At  (4.27a) 

(  iM(t)  >  =  Pj(t)  =  P,(0)  c'At  k4.27b) 


where  Pn  is  the  probability  of  being  in  state  n,  and  it  is  assumed  that  both  excited 
states  have  a  spontaneous  emission  rate  A.  The  middle  two  cases  are  equivalent, 
representing  complex  conjugates  of  each  other.  Taking  the  expectation  value  of 
equation  (4.21), 


(  ^12(0  )  —  P21  e  e 


(4.28) 


Note,  this  is  the  same  for  both  the  linear  and  quadratic  cases;  only  the  expressions 
for  a>12(t)  is  different. 

Finally,  using  equations  (4.27)  and  (4.28),  the  equation  for  power  radiated 
becomes 


P  = 


44oJff2QL2p2  +  l4oJgiaJ2pt  + 


3c* 


*  H  */a,w(t')dt' 
/h0'/i20^2le  ° 


+  fho-PioPitf  0 


-At 


(4.29) 


Notice  this  is  very  similar  to  equation  (2.21)  in  the  original  description  of  quantum 
beats,  except  for  the  exponential  terms  which  are  now  dependent  on  the 
instantaneous  transition  frequency,  w12(t).  This  reflects  the  Townes-Merritt  effect 
in  the  atomic  energy  states,  and  will  give  rise  to  beat  signals  that  depend  on  the  RF 
field.  Again,  the  distinction  between  the  linear  and  quadratic  cases  is  the  expression 
for  a/12(t). 
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Linear  Stark  Shift.  For  a  linear  Stark  shift  of  a  single  excited  level,  the 
instantaneous  energy  difference  is  represented  by  the  frequency 


W12(t)  =  W12  —  flyCOS  iA 


(4.30) 


The  exponential  term  in  equation  (4.29)  then  becomes 


i f  a>i2(t')  dt' 
e 


=  e 


iu>12t 


— i-h'sin  iA 
e 


(4.31) 


The  last  exponential  term  in  equation  (4.31)  can  be  expanded  into  a  Bessel  function 
series: 


if  w12(t')dt' 
e  0 


i{ul2— ni/)t 
e 


(4.32) 


The  expression  for  radiated  power,  equation  (4.29),  is  then 


P  = 


4^ig”|3pa + J>(a,) 


n ,  *(^12— n£/)t 


*  H  ^  O  *'(W12 

+  'M'^)  e 


n  =  'Oo 


-At 


(4.33) 


The  series  of  exponential  terms  indicate  the  possibility  of  quantum  beats  at 
frequencies  ul2  —  n v  depending  on  the  magnitude  of  the  coefficients.  To  illustrate 
this,  assume  a>10  a  u20  a  a>0  and  fto  =  fho  =  ^0-  addition,  let  the  initial 

H  H 

excitation  result  in  P,  =  P2  =  p{2  =  p2{  .  Then, 
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P  =  2Poe  At  [l  +  £:„(§")  cos(Ul2-n.-)t]  (4.34) 

n 

4^4  I  I  2 

where  P0  =  — §^3  •  Expanding  the  series  to  a  few  terms  around  n  =  0, 

P  8  ‘2P 0e  At[l  +  J0(§»)cosw12t 

+  Ji(p‘,)a’s(wI2-i/)t  +  •I-r(§v)cos(“i2+1')t  ■  •  ■]  H-33) 

This  shows  the  presence  of  beats  with  frequencies  u12  ,  ul2  ±  v  , .  . . ,  with  relative 
magnitude  J,,^)  •  Since  one  expects  v  >  Clv  ,  only  a  few  terms  around  n  =  0  will 
be  significant,  and  the  greatest  contribution  is  likely  to  come  from  the  J0  term, 
which  corresponds  to  the  beat  frequency  of  the  unperturbed  atom.  To  detect  the 
other  beat  frequencies,  consider  experimental  apparatus.  In  measuring  the  decay 
curve,  the  equipment  will  be  limited  by  an  effective  bandwidth,  BW.  Only  those 
signals  which  vary  at  a  frequency  less  than  this  bandwidth  can  be  measured.  Any 
beat  frequency  which  satisfies  w12  +  nv  >  BW  ,  will  average  to,  resulting  in  no 
impact  on  the  decay  measurement.  If,  however,  for  some  value  m,  u]2  -  mv  <  BW  , 
the  beat  signal  would  be  detectable  and  would  have  magnitude  J^-^)  .  The 
radiative  decay  would  then  follow 

P  =  2P„e~At  (1  +  Jm(9")  cos(u>]2— rm/)t)  (4.36) 

Another  possible  means  exists  for  a  very  large  field  such  that  the  term  J0(^)  «  0  . 

In  this  case,  the  higher  order  Bessel  functions  would  be  significant,  representing 
various  beat  signals.  These  experimental  considerations  will  be  covered  in  more 
detail  later. 
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Dressed  States  Calculation.  The  linear  Stark  shift  provides  a  relatively 
simple  case  which  can  be  examined  with  the  dressed  states  model,  without  great 
mathematical  complexity.  This  allows  another  comparison  to  the  adiabatic 
approach,  and  further  validates  that  model.  For  the  dressed  states  calculation, 
consider  the  energy  level  diagram  of  figure  13.  Beat  signals  can  occur  when  there  is 
the  possibility  of  transitions  from  a  level  in  each  of  the  two  excited  state  ladders  to 
a  common  level  in  the  ground  state  ladder.  In  figure  13,  the  j  l,n  )  to  |0,m  ) 
transition  could  interfere  with  the  |2,p  >  to  |0,m  )  transition  producing  a  beat 
signal.  Of  course,  in  the  dressed  state  picture,  this  combination  is  summed  over  the 
entire  photon  ladder,  so  the  expression  for  the  radiated  power  from  equation  (2.20), 

P=3^y  (4.37) 

n  ^p>m 

is  summed  over  the  combined  atom— field  states.  The  n  and  p  indices  each  refer  to 
the  excited  state  ladders  and  m  represents  the  ground  state  ladder. 

When  performing  the  summation,  there  are  two  cases  of  interest  for  the 
dipole  operator,  transitions  from  the  first  excited  atomic  state,  which  is  Stark 
shifted,  to  the  ground  state,  and  those  from  the  second  excited  state  to  the  ground 
state.  The  dipole  matrices  can  be  simplified,  using  equation  (3.110)  to  define  the 
energy  states,  resulting  in 

<“2ji|/ij0,m  )  =  fao <5nm  (4.38a) 

<“TjT|(i|0,m  >  =  MioJm-n(^)  (4-38b) 


With  the  approximation  u;nm  =  <jpm  =  ui0  ,  equation  (4.37)  becomes 
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p  =  ^^Ji^«|\^„pUei(,1'p)'/t+  ift.i2v„vP4.^i(n“p,'/t 

"'Vft0.^owpin,!!pl„eila'12  +  (n'pWt 

+  ho-for.tr.  +  <n'P)l',‘}  ^  (4-39) 

where  the  term  /?2pln  represent  the  density  matrix  element  between  states  |2,p  ) 
and  |  l,n  ).  The  first  term  in  the  summation  reduces  to 

^l/*2ol  P2m2m=  ^  I /*20 1  =  I  020  I  ^2  (4.40) 

m  m 

where  is  the  probability  of  being  in  state  |2,m  ).  Using  the  relation 
£  Jm-nJm-D  =  ^nD  [9]  produces  a  similar  result  in  the  second  term  of  the  summation: 

^l/ho!  Plnln  ~  Iftol  ^1  (^-41) 

n 

This  leaves  the  last  two  terms,  which  arc  similar  to  the  last  two  terms  of  equation 
(4.33).  The  primary  difference  is  in  the  density  matrix  /?2pln  which  depends  on  the 
probability  distribution  in  the  two  excited  state  ladders.  To  show  equivalence 
between  equations  (4.33)  and  (4.39),  first  change  the  summation  indices  such  that 
the  combinations  are  summed  over  all  decay  pairs  q  with  a  constant  difference  in 
photon  state  r  =  p  —  n: 

n  ,  P 

£2>o-/4v5„.  r,/)t 

r=  -00  a 


(4.42) 


Consider  now  the  density  matrix  summation,  specifically  the  term  Jr  E  P2q*riq  • 
Since  the  Bessel  fun  :tion  argument  is  small,  the  higher  order  Bessel  functions,  |  r  |  > 
1,  are  negligible.  Thus,  the  summation  over  q  need  only  consider  a  limited  range  of 
r.  For  a  typical  photon  distribution  in  a  real  field,  the  coefficients  of  nearby  number 
states  will  be  approximately  equal,  so 


Jr^/^qniq  *  *4^02. 


P2qlq 


(4.43) 


If 

This  last  sum  is  the  equivalent  of  p21 .  Equation  (4.39)  then  becomes 


P  ~  3c^''  ’tool  p2  +  l/hol  pi  +  /ho'/*20  02i^  e^Wj2  TU^ 


n  =  -oo 


+  *•**. 


r- -oo 


—At 


(4.44) 


which  is  of  the  same  form  as  the  adiabatic  result,  equation  (4.33). 

The  dressed  states  calculation  agrees  with  the  adiabatic  approach  when  an 
assumption  is  made  concerning  the  photon  state  distribution.  The  assumption  of 
equal  populations  in  nearby  states  is  quite  reasonable  for  real  fields.  This  agreement 
gives  greater  confidence  in  the  adiabatic  approach,  which  is  far  simpler,  particularly 
for  a  quadratic  shift  or  for  level  crossing  calculations. 


Quadratic  Stark  Shift.  When  a  quadratic  Stark  shift  is  present,  the 
instantaneous  energy  difference  between  upper  levels  is  represented  by  the  frequency 
given  in  equation  (4.22) 


o>12(t)  =  u)n  4-  0'cos2</t 


(4.45) 
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where  ft'  =  U^1  -10 .  Notice  that  if  ft.  «  ft,  ,  then  ft'  «  0  and  there  would  be 

CJi  1  * 

no  beating  effect.  Similar  to  the  development  in  the  linear  case,  the  exponential 
term  in  equation  (4.30)  becomes 


dt'= 


=  e 


O'  O'  . 

“b  h^- 1  ■+■  i^^isin  2iA 


=  M 


O' 

q  i  i(ojj2+  ^  +2nt/)t 


(4.46) 


This  is  similar  to  equation  (4.32)  in  the  linear  case,  except  the  beat  frequencies  are 

o' 

at  even  harmonics  of  i/,  given  by  ui2  +  +  2nv .  Also,  there  is  a  slight  shift  in 

O' 

the  fundamental  frequency  by  the  amount  y . 

Using  equation  (4.46),  the  power  radiated  is 


P  = 


4A|ffo|2p,  +  ^Aigiolip, 


3? 


,  4uMoL  *  .«Vt  /fl\ 
+  ^20  P 21 


O' 


Qi  t'(u>l2+  Y  +2ni/)t 


+  th O'PlO  ^2^Jn(^)  e 


ft' 


,(a’i2+  y  +2m/)t-|1 


—At 


(4.47) 


Beats  at  frequencies  other  than  w,2  are  expected  to  be  much  smaller  than  the  linear 
case  since  ft'  is  small  compared  to  ft,,.  The  dominant  term  in  the  expansion  is  then 

ft' 

J°("21/)  '  eva^uate  the  possibility  of  beats  at  other  frequencies,  first  consider  ft'. 
For  ft'  to  be  non— zero,  the  individual  Rabi  frequencies  ft,  and  ft2  can  not  be  equal, 
so  /iw  t  fro  .  In  this  case,  assume  the  dipole  moments  differ  by  a  constant,  or 
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H  H 

P2 o  =  K0,o.  Also,  as  before,  let  ul0  a  and  Pj  =  P2  =  pn  =  P21  •  Then, 
equation  (4.47)  becomes 


P  = 


P0e 


-At 


1  +  K2  -f-  2K^T  Jn(^)  cos(wj2  +  ^  +  2nt/)tl  (4.48) 


To  observe  beats  resulting  from  the  Townes— Merritt  interaction,  consider 
again  the  effective  bandwidth  BW  .  If  there  is  a  value  m  for  which 


2  4*  ^  •+•  2mi/  <  BW 


(4.49)  • 


then  an  amplitude  modulation  in  the  decay  curve  at  that  frequency  could  be 
observed.  The  power  radiated  would  then  be 


P  =  Poe  At  [l  +  K2  +  2KJm(^)  cos(u>j2  +  S'  +  2ni/)t, 
► 


(4.50) 


It  is  anticipated  that  this  modulation  would  be  much  less  than  in  linear  case 
due  to  a  small  IP  and  the  factor  K.  Here  the  modulation  amplitude  would  be 


2K  -Ufj) 

1  +  R'2  “ 


(4.51) 


as  opposed  to  the  linear  case  where 

Ml  =  J„(9")  (4  W) 


Nonetheless,  equations  (4.36)  and  (4.50)  demonstrate  the  possibility  of  quantum 
beats  which  result  from  interaction  with  the  RF  field. 
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The  additional  beat  signals  result  from  an  effective  splitting  of  the  excited 
states  due  to  the  nonresonant  RF  field.  When  any  one  of  the  levels  in  an  excited 
state  ladder  comes  close  to  a  level  in  the  other  excited  state  ladder,  a  beat  signal  in 
the  radiative  decay  could  be  observed.  The  magnitude  of  this  modulation  depends 
primarily  on  the  RF  field  strength  and  the  experimental  setup,  as  will  be  discussed 
in  the  next  chapter. 

Level  Crossings 

If  the  two  excited  states  have  the  same  energy,  a  level  crossing  condition 
exists,  whether  the  crossing  is  the  normal  occurrence  or  a  result  of  an  energy  shift  of 
one  or  both  levels.  When  a  crossing  occurs,  the  two  possible  transitions  can 
interfere  causing  a  shift  in  the  atomic  fluorescence  pattern,  as  compared  to  the 
uncrossed  situation.  The  interference  is  heavily  dependent  on  the  specific  transitions 
involved  and  the  observation  angle.  From  equation  (2.23),  with  the  atom  at  the 
origin,  the  radiated  energy  as  a  function  of  the  point  of  observation  is  given  by 

s = (r1i&iX|1  X  <W‘>  <4-53> 

n  ^>i» 

Note  the  Einstein  summation  notation  is  used  for  the  i  and  j  indices.  This  section 
will  explore  level  crossings  resulting  from  an  interaction  with  the  RF  field.  As 
previously  discussed,  the  RF  field  causes  an  effective  splitting  of  the  atomic  levels, 
which  is  easily  seen  with  the  dressed  states  model.  However,  for  mathematical 
simplicity,  the  calculations  here  will  follow  the  adiabatic  approach  in  defining  the 
atomic  dynamics,  which  has  compared  favorably  with  the  dressed  states  result  in 
simpler  cases.  Both  the  linear  and  the  quadratic  Stark  shift  will  be  covered,  with 
equations  (4.12)  and  (4.21)  used  to  describe  the  atom  in  the  RF  field  for  their 
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respective  cases.  The  approach  is  similar  to  that  used  in  the  case  of  quantum  beats. 
For  the  three  level  atom,  the  sum  in  equation  (4.53)  reduces  to  four  terms: 


Ao  /4op2e  "h  /%/4oPie 
+  Ao  A0P2M  +  Ao  AoPM 


(4.54) 


assuming  tu10  «  u>2() .  In  evaluating  this  equation,  the  geometry  of  the  experiment 
and  the  dipole  moments  of  the  specific  transitions  involved  become  very  important. 
With  regard  to  the  experimental  geometry,  the  orientation  of  the  excitation  pulse 
polarization  and  the  viewing  angle  are  critical  parameters.  Geometries  which  yield 
the  greatest  effects  will  be  considered  here.  The  transition  dipole  moments  impact 
the  degree  of  interference,  in  much  the  same  way  as  the  orientation  of  radiating 
elements  in  an  array  determine  the  overall  antenna  pattern.  Three  possible 
combinations  for  the  two  transitions  will  be  covered. 

For  a  dipole  transition,  the  moment  can  be  written  as 


(  n  |plm  )  =/W 


(4.55) 


where  €  is  a  unit  vector  in  the  direction  of  /*.  The  orientation  of  e  influences  the 
initial  excitation  and  determines  the  polarization  of  the  emitted  radiation.  This 
orientation  depends  on  the  initial  and  final  orbital  quantum  numbers,  1,  and 
magnetic  quantum  numbers,  m  [8].  Table  1  summarizes  the  possibilities.  For  the 
transition  in  which  the  magnetic  quantum  numbers  remain  the  same,  the  dipole 
moment  is  oriented  in  the  z  direction  and  radiation  would  then  be  linearly  polarized 
in  the  same  direction  if  view  in  the  xy  plane.  The  other  cases  result  in  circular 
polarization  if  viewed  along  the  z  axis.  Note  that  even  though  these  transitions  are 
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Quantum  Numbers 


Orbital 

Magnetic 

Dipole  Moment 

r  =  1  ±  1 

m1  =  m 

r  =  U1 

m1  =  m  —  1 

^nm(x  +  *y)/^  =  /W+ 

r  =  u  i 

m1  =  m  +  1 

/'na/*  -  *)/&  E  /W- 

Table  1.  Transition  Dipole  Moments.  Initial  state  is  defined  by  1  and  m. 

considered  orthogonal,  at  other  viewing  angles,  the  radiation  from  these  transitions 
do  not  have  a  simple  polarization,  and  thus  may  interfere  with  each  other.  With 
regard  to  the  initial  excitation,  the  orientation  of  the  excitation  pulse,  relative  to 
the  dipole  moments,  determines  the  initial  values  of  the  density  matrix  elements. 

For  example,  to  excite  a  transition  with  the  dipole  moment  in  the  z  direction,  the 
pulse  should  optimally  be  polarized  in  the  direction. 

When  the  atom  decays,  the  degree  of  interference  depends  on  the  -elative 
orientation  of  the  dipole  moments.  The  three  case3  tr-  be  considered  are 

1)  .  t. 

2)  c,  ,  z 

3)  z  ,  z 

These  are  representative  of  all  possible  combinations,  and  ea^n  case  will  be 
considered  for  both  the  linear  and  the  quadratic  sh::  s.  In  each  of  these  case,  an 
optimal  orient?*  ion  for  the  excitation  pulse  and  the  viewing  angle  will  be  specified. 

Linear  Stark  Shift  As  in  the  case  of  quantum  beats,  the  diagonal  elements 
of  the  density  matrix  represent  the  probability  of  the  atom  being  in  ihe  given  energy 
state.  From  equation  ( '.27),  these  are 


!•) 


(4.56a) 

(4.56b) 


P22M  —  ^2e 

Pii(t)  =  P,e“At 


The  off  diagonal  elements  p2l  and  p12  represent  the  interference,  and  are  critical 
to  the  level  crossing  effect.  These  are  complex  conjugates,  so  only  one  will  be 
treated  explicitly.  Recall  from  equation  (4.28), 


Pull)  =  Pi i  e  Ate 


(4.57) 


For  a  linear  Stark  shift,  the  instantaneous  frequency  difference  is 

w12(t)  =  u>12  —  (lvcos  Vi  .  Using  equation  (4.32),  the  second  exponential  term  then 

reduces  to  a  Bessel  function  series: 


~  ^21  e  At  Jn(^y)  e 


n  i(un-nv)t 


(4.58) 


Equation  (4.54)  then  becomes 


,n ,  '(w.2-n")t 


S  =  (rZfe»HXl)^Uo  PioP*  +  A*  /*!>,  +  Am  e 

n 

.*  H  V1  O  ~ t(^i2— 1 nJ/)ti  a  * 

+  fib  e  e  (4  59) 


Further  ev equation  requires  expressions  for  the  the  dipole  moments  and  the  initial 
density  matrix,  whic.  :epend  on  the  specific  combination  of  dipole  moments.  The 
three  cases  previously  discussed  will  be  covered  here. 
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Case  1:  First  let  the  two  dipole  moments  be  given  by 

fh  O  =  /x0e+  =  ^o(x  +  v)/&  (4.60a) 

tho  =  t*0e-  =  /<0(x  “  *y)/|2  (4.60b) 

Also,  define  the  position  coordinates  using  the  standard  spherical  coordinate  system: 

Xj  =  r  sin  0  cos  <(> 

x2  =  r  sin  0  sin  <f>  (4-61)  • 

x3  =  r  cos  0 

With  these  expressions,  the  summations  over  i  and  j  in  equation  (4.59)  can  be 
reduced.  Table  2  lists  the  individual  terms  involved  in  the  summations. 

VWo  =  6i)4o4o  =  & 

.*  .* 

^ij/h'o^O  =  =  0 

Xi^^io = sin7° 
xixj/4o/4o  =T-^s\n20 
xixj/4/4o  =  ^sin^e12^ 

xixj/4o/4o  =  T-^  sin2^ 

Table  2.  Dipole  Moment  Products  for  Case  1. 

The  radiated  power  as  a  function  of  observation  angle  is  then 
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.,2 .  .4 
Q  _  fipUfo 

4?rc3r4 


(1  +  cos2«)(P2  +  P.)  -sin2«  [c^^J^) 


+  e 


n 


n„x  e-*(^i2-ni/)t 


-At 


i(a;12— n^)t 
e 


(4.62) 


Next,  expressions  for  the  initial  density  matrix  elements  are  required.  These  depend 
on  the  orientation  of  the  excitation  pulse.  An  element  of  the  density  matrix  can  be 
written  as 

4  =  Vp  (4-63) 


where  an  is  the  wave  function  coefficient  corresponding  to  the  atomic  state  |  n  ) 
after  initial  excitation.  For  the  excited  states,  an  is  given  by 


tC'fi noQ 


(4.64) 


where  e  is  the  polarization  of  the  excitation  pulse,  and  Q  is  the  area  of  the  pulse, 
specifically  given  by  Q  =  ^-j£(t)dt  .  This  assumes  a  short  pulse  of  frequency  u  and 
duration  T,  such  that  (un0  -  w)T  <  1  . 

Let  the  excitation  pulse  be  linearly  polarized  in  the  y  direction,  so  e  =  y, 
which  will  equally  excite  both  transitions.  Figure  14  shows  the  orientation  for  this 
case.  The  initial  values  for  the  density  matrix  are  then 


2  2 

Pj  =  P2^^S- 
2  2 

H  H  ll0  Q 

Pi  1  =  P 12  =  ~iTL- 


(4.65a) 

(4.65b) 
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Figure  14.  Level  Crossing  Geometry 
Substituting  the  density  matrix  terms  into  equation  (4.62)  yields 


S  =  e  At  (l+cos20)  +  sin20^Jn(^)  cos(u;12t-oi/t-20)  (4.66) 


n 

Integrating  S  over  time  gives  the  total  energy  radiated  by  the  atom  in  the 
direction  specified  by  0  and  (f>: 


J  =  f*S(t)  dt 


/q/41 
!  4  jrc3  r 


1+C0S2# 


+s\n2ey  Jn(^)Acos  2^_-f(^2-nWsin  2 f 

A  Zu  n  ^  (u>,  2-n^r+A2 


(4-67) 


This  result  is  useful  in  predicting  the  radiated  energy  as  a  function  of  the  RF 
frequency,  u  .  Level  crossings  occur  when  ul2  -  nv  =  0  ,  which  agrees  with  the 
dressed  states  picture.  When  ?.  level  crossing  condition  exists,  there  is  a  maximum 
in  the  interference  term,  depending  on  the  observation  angle,  with  relative 
magnitude  Jn(^)  . 
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4  3  2 

Figure  15.  Level  Crossing  Signal  for  a  Linear  Stark  Shift 


For  an  example  of  a  possible  experimental  geometry,  let  the  observation 
angle  be  defined  by  6  =  | ,  <f>  =  0  .  Then 


J~J°  1  +  Z  TlST^TP1 


(4.68) 


Whete  Jo  5  and  An  =  ((J12— m/)2 .  This  indicates  a  series  of  peaks  in  the 

radiated  energy  as  v  changes.  Figure  15  shows  a  representative  level  crossing  signal 


for  this  case.  When  mv  «  ,  then  for  n  t  m  ,  the  interference  term  is  negligible 

since  ^  <  1  : 


1  l  J  A2 

1  +  KTTJP 


(4.69) 


This  represents  the  signal  strength  near  one  of  the  peaks,  and  is  similar  to  the 
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X 


original  level  crossing  calculation,  equation  (2.29).  The  magnitude  of  the  peak  is 
Jmax  =  J0[  1  +  Jm(^1/)]  and  the  width  is  approximately  ^  .  This  becomes  the  basis 
for  using  a  nonresonant  RF  field  as  a  spectroscopic  tool.  If  the  peaks  can  be 
detected,  then  the  values  of  ^at  which  the  crossings  occur  determine  the  original 
energy  level  spacing.  Note  that  other  viewing  angles  will  produce  different  types  of 
level  crossing  curves  than  figure  15. 

Case  2:  When  the  relevant  transitions  have  dipole  moments  given  by 

/ho  =  »oc+  (4.70a) 

fho  =  /V*  (4.70b) 

the  potential  interference  effect  is  expected  to  be  diminished.  This  compares  to 
antennae  placed  along  orthogonal  axes;  there  is  limited  overlap  between  the  beam 
patterns.  Table  3  lists  the  dipole  products  for  this  case. 

•  *  .* 

.  *  $ 

V!o/4o =  = 0 

XiVlo/4  =  sin20 
.  * 

xixj/4/4o  =  rVo  cos20 

xix  /iJo/4o  =  ^ ^  sin  2 9e1^ 

2 (2 

XjXj/t^oHo  =  1 —  sin  20  e 
2 |2 

Table  3.  Dipole  Moment  Products  for  Case  2. 

Substitution  of  the  di j>olo  products  into  equation  (4.59)  results  in 
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s  —  jt 2yo. 

b  ”  2xc3r4 


B-'(l  -co^P,  +  (1  _sm!i)Pi_?inW 


2(2 


n 


p'VJ.X-I.tv-)  e 
] 


f)  *(W|2-n/y)t 


(4.71) 


This  shows  that  there  is  no  interference  for  9  =  C,^  ,  which  is  consistent  with  the 
antenna  analogy.  Further,  to  obtain  excitation  of  both  levels,  the  initial  pulse 
would  have  to  be  linearly  polarized  at  an  angle  between  the  z  axis  and  the  xy  plane. 

For  the  excitation  pulse,  let  e  =  (y  +  z)/(S  ,  a  45  degree  angle  to  the  z  axis. 
The  initial  values  for  the  density  matrix  are  then 


P 


i  ~ 


M&2 


(4.7‘2a) 


p  _  uoQ 2 
r2  -  2^ 


4£2 

2{T 


4a2 

2(T 


(4.72b) 

(4.72c) 

(4.72d) 


The  expression  for  radiated  power  then  becomes 


e  At  (4.73) 


Integration  results  in 


o  _  /4<4Q2 

4xc^r  4 


1  ,  3sin20  sin  29 
K  *r  7 - 7 


^Jn(^)  sin(w12t-ni4+^) 


J  = 


^o^toQ2 
'  xc3  r 1 


2-f3s  in  29 


IX 


sin20 


(4.74) 


which  is  similar  to  equation  (4.67)  in  case  1.  Again,  maxima  will  occur  whenever 
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ul2  -  nv  =  0  ,  coinciding  with  the  level  crossings. 

If  the  detector  is  placed  at  0  =  j  ,  <j>  =—  ^  ,  the  total  energy  received  is 
given  by 


3 -30 


7,1 

S  +  ? 


A2 

A2  +  A'i 


(4.75) 


The  result  compares  to  the  previous  case,  except  that  the  interference  effect  is 
reduced  by  approximately  ^ .  The  energy  received  as  a  function  of  v  would  look 
like  figure  15,  but  the  magnitude  of  the  peaks  is  given  by 


'max 


(4.76) 


The  unfavorable  geometry  is  the  reason  for  the  lesser  effect. 


Case  3:  The  most  favorable  geometry  exists  when  the  transition 
dipoles  have  the  same  orientation: 


Mio  =  /^20  =  Voz  (4."~ ' 

This  results  in  a  maximum  overJap  in  the  radiation  patterns.  Using  the  dipole 
products  given  in  table  4,  the  expression  for  radiated  energy  becomes 

S  = 


Notice  there  is  no  <fi  dependence. 


al<4>  „-At 


+  e 


P*  +  P.  +  e 

n 

n  -t(w12-ni/)t 


0  i(w12- n^)t 

e 

sin2# 


(4-71) 
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$  ★ 

VMo  =  ^ij^oMo  =  & 

*ijMio/4o  =  SijMd  =  Mo 

xixjMloMio  =  xixj/4oMo  =  f2Mo  cos2^ 


.* 


=  xixjM2oMio  =  rV2  cos2# 


Table  4.  Dipole  Moment  Products  for  Case  3. 
The  optimum  excitation  polarization  is  c  =  z  ,  in  which  case 


P,  —  P2  —  p<i\  —  P12  —  M2Q2 


(4.79) 


Equation  (4.78)  then  becomes 


S  = 


1  + 


£  J^)  cos(cu12-m/)t 


sin2# 


(4.80) 


Once  again,  integration  results  in 


J  = 


-4-4o21 

xc^r  • 


IJn(^)A-2  T iq 


sin2# 


(4.81) 


It  is  obvious  that  the  maximum  radiation  occurs  for  0  —  ^ .  As  in  the 
previous  two  cases,  peaks  occur  when  An  =  w12  —  nu  =  0  ,  and  in  this  case  have 
magnitude 

=  OJl  +  -U^)]  f4-82) 


This  compares  favorably  with  the  first  case,  having  the  same  relative  magnitudes 
between  the  peaks  and  the  background,  but  this  case  has  a  larger  absolute  value. 
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Table  5.  Level  Crossing  Signals  for  a  Linear  Stark  Shift 
All  three  cases  for  an  atom  with  a  linear  Stark  shift  of  an  excited  state 


exhibit  changes  in  the  fluorescence  resulting  from  the  RF  field.  Table  5  summarizes 
the  results  for  each  case  with  a  specific  geometry.  Recall,  J0  h  ^4^3^  ~  and  An  = 
ul2  —  nv.  In  these  results,  the  presence  of  a  maximum  indicates  a  level  crossing 
condition.  Other  observation  angles  can  result  in  different  level  crossing  signals,  for 
example  a  decease  in  the  energy  received  as  the  levels  cross.  In  any  case,  if  these 
can  be  measured,  then  the  original  energy  separation  can  be  determined  by  the  RF 
frequencies  at  which  crossings  occurred.  The  magnitude  of  the  Bessel  function 
terms,  Jn(^)  ,  will  be  key  in  determining  whether  the  peaks  can  be  detected. 
Overall  though,  the  linear  case  presents  the  best  opportunity  for  observing  the  level 
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crossings,  when  compared  to  quadratic  shifts,  as  will  be  shown. 

Quadratic  Stark  Shift.  For  an  atom  with  a  quadratic  Stark  shift,  the  off 
diagonal  density  matrix  elements  are  given  by 


H  -At 

P21W  —  P21  e  e 


(4.83) 


where  u>12(t)  =  u;,2  +  Q'cos 2iA  and  O'  =  t^ie  Bessel  function 

expansion,  this  becomes 


/ ,  \  H  — A.  t 

P2M  “  P 21  e 


X,  ,n\  ‘("u+  ?  +2ni/)t 
Jn<5P> e 
n 


(4.84) 


Notice  that  this  is  of  tiie  same  form  as  the  linear  Stark  case,  equation  (4.58),  except 
for  the  Bessel  function  argument  and  that  only  even  harmonics  of  the  RF  frequency 
are  included.  As  in  the  section  on  quantum  beats,  if  fl,  ks  fl2  ,  then  fi'  s  0  and 
there  would  be  a  negligible  effect  at  level  crossing,  since  the  only  significant  term 

O' 

would  be  J 0(-^)  *  1  .  Thus,  once  again,  it  is  assumed  the  magnitudes  of  the  dipole 
moments  differ  by  a  constant,  such  that  =  K  | /z10 (  . 

Substitution  of  equation  (4.84)  into  equation  (4.54)  yields  an  expression  for 
the  radiated  power: 


*A't 


s  =  /4*P2  +  /4’p,  +  /4/V2,XJ"C>  e  " 

n 

.  *  ii  t-i  o'  ~  *"A  nt"J  A  * 

+  ^20  tl\oP  12^^ e  C  (4.85) 
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I  Mio  =  /V* 
fho  =  /h>c- 


c  =  y 


J  =  ^[(l+cos20)(l+K4)+K2sin20]T  J  n(^)A^§2|±AA^sin2^ 

L  n  n  J 

0  =  ff/2  ,  <j)  =  0 


J  =  l 


o  1  +  R4  +  K2XJn(fe  A2  T  A;2. 


fto  =  /*0C+ 
tho  =  /xoz 

^  =  ^  f2+(4K4-l)sin2^ 


f  =  (y+z)/^ 


t 

K2sin2^V*  t  /fi\A2sin0+AA nc.os  6 
- 3 - 2/n(^} - A^T^T2 

n 


0  =  7r/4  ,  (p  —  — tt/2 


y  = 


f  f2+(4K*-l) 

*/o  8 


S’* 


A2 


III 


MlO  =  /*oz 

fho  =  W 


J  =  VC 


1  +  K4  +  2K2^Jn(^)  £7- 


+ 


IT2 

nn  J 


sin20 


0=  tt/2 


J  ~  2  Jr 


1  +  K4  +  2K*j;j.(j&  s 


A2 


Table  6.  Level  Crossing  Signals  for  a  Quadratic  Stark  Shift 

O' 

where  A,!,  =  u/12  +  +  2n^ .  This  is  comparable  to  equation  (4.59)  of  the  linear 

case;  hence,  the  results  for  the  three  combinations  previously  considered  would  be 
similar  here.  These  are  summarized  in  table  G. 
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Figure  16.  Level  Crossing  Signal  for  a  Quadratic  Stark  Shift 


In  each  case,  the  energy  received  varies  as  a  function  of  v,  as  a  result  of  level 

O' 

crossings.  Peaks  occur  when  A„  =  ul2  +  +  2nt/  =  0  ,  with  relative  maxima 

o  • 

determined  by  Jn('27;)  •  Figure  16  shows  the  received  energy  versus  RF  frequency 

for  case  1.  This  is  very  similar  to  the  linear  case,  except  the  peaks  occur  at  half  the 

frequency,  reflecting  the  even  harmonics  of  the  RF  frequency.  In  addition,  each 

O' 

peak  is  shifted  slightly  from  even  fractions  of  u2x  by  the  amount  ^  due  to  the  static 

A 

Stark  shift.  The  width  of  each  peak  is  — ,  again  making  this  potentially  useful  for 
high  resolution  spectroscopy.  A  noteworthy  point  is  that  this  technique  would  be 
relatively  doppler  free.  Normally,  the  doppler  effect  results  in  a  significant  line 
broadening;  however,  in  this  case,  the  level  shifts  result  from  interaction  with  the 
RF  field.  Doppler  shifts  at  RF  frequencies  are  very  small  compared  to  the 
transition  frequency. 

Limiting  factors  are  the  relative  magnitudes  of  the  dipole  moments,  indicated 
by  the  constant  K,  and  the  magnitude  of  the  Bessel  function  coefficients  Jn(^)  ■ 
These  will  determine  whether  the  peaks  can  be  detected,  and  the  technique  usable. 
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The  previous  two  chapters  presented  theoretical  descriptions  of  the 
Townes-Merritt  effect,  quantum  beats,  and  level  crossings,  which  was  the  goal  of 
this  dissertation.  As  a  final  step,  a  brief  look  at  some  experimental  aspects  will 
provide  further  insight  to  these  phenomena.  An  example  is  the  field  strengths 
required  to  observe  these  effects  for  a  given  atomic  sample.  The  intent  here  is  not 
to  describe  an  experiment,  but  to  give  an  idea  of  the  feasibility  of  such  an 
experiment  by  considering  order  of  magnitude  estimates  of  critical  parameters.  In 
the  Townes-Merritt  effect,  the  primary  concern  is  the  relative  magnitude  of  the 
sidebands.  Experimental  requirements  for  observing  the  sidebands,  in  linear  and 
quadratic  cases,  will  be  examined,  and  these  conditions  will  also  apply  to  level 
crossings  and  quantum  beats.  An  additional  concern  in  quantum  beats  is  the  beat 
frequency,  and  for  level  crossings,  the  viewing  angle  is  important. 

The  calculations  on  the  Townes-Merritt  effect  describe  absorption  sidebands 
at  integer  multiples  of  the  RF  frequency  for  a  linear  Stark  shift  and  even  integer 
multiples  for  a  quadratic  shift,  with  relative  amplitudes  of  J|!(— ")  and  > 

respectively.  In  order  to  observe  these  sidebands  there  must  be  sufficient  separation 
between  adjacent  signals  so  that  the  natural  linewidth  does  not  obscure  individual 
sidebands,  and  the  magnitudes  governed  by  the  Bessel  function  terms  must  be  large 
enough  for  detection.  A  condition  that  v  >  A  ensures  separation  of  individual 
sidebands.  Typical  values  for  A  are  on  the  order  of  106  — ♦  109  sec'1,  so  a  nominal 
value  for  v  is  10  times  A  or  107  — *  1010  Hz. 

The  amount  of  energy  in  the  sidebands  is  determined  by  the  Bessel  function 
coefficients.  For  a  given  sideband  to  be  observable,  the  magnitude  of  the  associated 
coefficient  should  be  on  the  order  of  Jn  a  .01  — » .1  as  a  minimum.  This  then 
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places  a  bound  on  the  Bessel  function  arguments  and  subsequently  defines 
experimental  constraints. 

The  key  term,  from  an  experimental  perspective,  in  the  Bessel  function 
arguments  is  the  Rabi  frequency,  ^ .  The  dipole  moment  for  an  allowed 
transition  can  be  approximated  as  /i  a  ea0  ,  where  a0  is  the  Bohr  radius, 
approximately  5  *  10'9cm.  The  Rabi  frequency  can  then  be  written 

nv  =  1075u  (5.1) 

where  £v  is  the  magnitude  of  the  RF  field  in  V/cm.  For  £v  =  100  — •  1000  V/cm  , 
=  io9  — » io10  Hz. 

Now  consider  a  specific  example  where  u  =  109  Hz  and  w0  =  5*  1015  Hz  . 
The  relative  intensities  of  the  sidebands  are  given  by 

Linear:  =  Jj(10%)  (5.2a) 

Quadratic:  =  (5.2b) 

This  illustrates  the  significant  decrease  in  effect  for  the  quadratic  case,  as  compared 
to  a  linear  Stark  shift.  For  an  RF  field  with  an  amplitude  of  100  to  1000  V/cm  the 
argument  for  the  Bessel  function  in  the  linear  case  is  between  1  and  10,  resulting  in 
a  significant  amount  of  energy  in  the  sidebands.  Table  7  lists  relative  intensities  for 
the  lower  order  sidebands  in  this  case. 

For  the  quadratic  case,  with  £v  =  1000  V/cm  ,  the  sidebands  are  negligible. 
The  first  sideband  has  a  relative  intensity  of  J^(  10-5)  «  2.5  *  IO-11.  To  observe 
sidebands  for  a  quadratic  shift  would  require  very  large  field  strengths  in  this  case, 
or  alternately  consider  transitions  with  lower  frequencies  such  that  the  Bessel 
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z 

•*o(z) 

J?(z) 

Jl(z) 

m 

Js(z) 

Ji(z) 

J?(z) 

1 

.5855 

.1936 

.0132 

.0003 

2 

.0501 

.3326 

.1244 

.0166 

.0011 

3 

.0676 

.1149 

.2362 

.0955 

.0174 

.0018 

4 

.1577 

.0043 

.1325 

.1850 

.0790 

.0174 

.0024 

5 

.0315 

.1073 

.0021 

.1331 

.1530 

.0681 

.0171 

.0028 

6 

.0226 

.0765 

.0589 

.0131 

.1279 

.1311 

.0604 

.0167 

7 

.0900 

.0000 

.  .0908 

.0280 

.0249 

.1210 

.1150 

.0545 

8 

.0294 

.0550 

.0127 

.0847 

.0111 

.0345 

.1139 

.1027 

9 

.0081 

.0601 

.0209 

.0327 

.0704 

.0030 

.0417 

.1072 

10 

.0604 

.0018 

.0648 

.0034 

.0482 

.0547 

.0002 

.0469 

Table  7.  Sideband  Intensities  for  a  Linear  Stark  Shift 


function  arguments  are  near  unity. 

In  determining  a  feasible  range  for  the  Townes— Merritt  effect  consider  the 
relationships  between  u>0  ,  A  ,  and  v  .  As  previously  discussed,  the  RF  frequency 
should  be  greater  than  the  atomic  decay  rate,  on  the  order  of  v  «  10A  ,  to  break  out 
individual  sidebands.  Further,  the  A  coefficient  is  dependent  on  the  transition 
frequency  [23]: 

A  =  (5’3) 


in  SI  units.  Again  assuming  /i «  ea0  ,  this  can  be  simplified  to  Am3*  10'4°u^  . 
Thus,  for  a  given  transition  frequency,  the  optimum  RF  frequency  can  be 
determined  and  subsequently,  as  in  equations  (5.2),  the  required  field  strength. 

As  an  example,  consider  the  hydrogen  atom  IS  to  2P  transition.  The  2P 
level  exhibits  a  linear  Stark  shift  due  to  mixing  between  the  2S  and  2P  states 
[32],  resulting  in  energy  shifts  for  the  two  states  of  ±3ea0<f.  The  transition  frequency 
is  approximately  u0  «  1016Hz  and  the  decay  rate  is  about  6  *  108Hz  .  A 
reasonable  choice  for  the  RF  frequency  is  then  v  =  3  *  109Hz.  The  relative 
intensities  of  the  sidebands  in  this  case  are  given  by  equation  (5.2a),  J2(10"2£„). 
Table  7  lists  the  intensities  for  typical  values  of  £v. 
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Figure  17.  RF  Field  Strength  Required  for  a  Quadratic  Shift 

For  the  quadratic  case  the  problem  is  to  determine  the  approximate  range 
where  the  effect  could  be  observed  for  a  given  maximum  field  strength,  representing 
an  experimental  constraint.  Assuming  v  ~  10A  ,  the  Bessel  function  argument  can 
be  rewritten  using  equation  (5.3) 


Jll  _  3a~c0c3 
2lj0u  20  hu$ 


el 

S3 


(5.4) 


02 

Considering  only  the  first  sideband,  let  J  >  01  be  the  criteria  for 

detection.  For  small  arguments,  much  less  than  one,  the  Bessel  function  can  be 
approximated  as  [9] 


^n(z) 


0fe/21” 


(5.5) 


Then  if  E  —  1000  V/cm  ,  using  equations  (5.4)  and  (5.5),  the  transition  frequency 
must  be  u0  <  5  *  10,4Hz  .  Following  the  same  approach,  figure  17  shows  the 
minimum  required  field  strength  for  a  range  of  typical  transition  frequencies. 
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Higher  field  strengths  would  allow  higher  order  sidebands  to  be  observed.  Thus, 
based  on  these  order  of  magnitude  estimates,  the  effect  can  be  observed  in  a  wide 
range  of  potential  samples,  shown  here  to  have  transitions  from  about  the  middle  of 
the  visible  into  the  infrared. 

Another  case  of  interest  involves  an  atom  which  exhibits  a  Zeeman  shift  in 
the  energy  levels.  This  is  a  linear  effect  proportional  to  the  magnetic  field  strength 
B.  The  interaction  Hamiltonian  for  the  Zeeman  effect  is 

H'E-ivB  (5.6) 

where  ^  is  the  magnetic  moment  and  B  is  an  oscillating  magnetic  field. 
Development  of  the  equations  of  motion  is  analogous  to  the  linear  Stark  shift.  For  a 
two  level  atom  in  which  only  the  excited  state  is  shifted,  the  interaction 
Hamiltonian  can  be  written 


The  instantaneous  transition  frequency  is  then 

uj'(t)  =  u0  -  Dbcos  iA  (5.8) 

where  Db  =  .  As  in  the  case  of  the  linear  Stark  effect,  there  are  absorption 

sidebands  at  u0  +  w  with  relative  intensity  J„(— )  • 

To  estimate  the  strength  of  the  sidebands  in  this  case,  the  magnetic  dipole 
can  be  approximated  by  the  Bohr  magneton,  /im  «  ^  .  The  intensities  then 
depend  directly  on  the  B  field  strength: 


93 


=  JS(2S75)  “  J2(10“  ®)  (5-9) 

where  B  is  in  webers/m2.  The  primary  experimental  constraint  arises  in  producing 
a  strong  B  field  at  the  relatively  high  frequency  v.  The  coil  inductance  limits  the 
maximum  field  strength  at  a  given  frequency. 

To  analyze  this  constraint,  consider  a  typical  setup  in  which  an  AC  power 
supply  with  a  maximum  voltage  V  drives  a  wire  coil  to  produce  a  B  field  at 
frequency  v.  The  field  strength  is  given  by 

B  =  /i0IjSin  (5.10) 

where  n0  is  the  permeability  of  free  space,  I  the  instantaneous  current,  N  the 
number  of  turns,  /  the  length  of  the  coil,  and  0m  the  angle  from  the  center  to  the 
edge  of  the  coil  [33].  The  current  flow  is  limited  by  the  maximum  available  voltage 
V  and  the  coil  impedance  XL  according  to  Ohm's  Law:  I  =  V/XL  .  The  impedance 
is  equal  to  the  product  of  the  field  frequency  and  the  coil  inductance,  which  is 

L  =  njj  A  (5-11) 

where  A  is  the  cross-sectional  area  of  the  coil.  Thus  the  maximum  B  field  strength 
can  be  expressed  as 

B  =  HTO5in(,»  <512> 

For  a  power  supply  with  a  maximum  voltage  of  104V  ,  and  a  coil  which  has 
N  =  100  ,  A  =  100  cm2,  and  0m  =  30°,  the  sideband  intensities  are  given  by 
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=  J^) 


(5.13) 


In  order  for  the  sidebands  to  be  observable,  the  Bessel  function  argument  should 
again  be  of  order  unity,  so  the  maximum  B  field  frequency  is  v  ss  107  — *  108  Hz  . 
The  atomic  decay  coefficient  is  then  an  order  of  magnitude  lower,  A  «  106  — >  107 
Hz  ,  and  using  equation  (5.3),  the  maximum  transition  frequency  is  on  the  order  of 
10' 5  Hz  ,  which  is  comparable  to  the  quadratic  Stark  shift  case. 

These  order  of  magnitude  calculations  demonstrate  the  experimental 
constraints  on  observing  the  Townes— Merritt  effect  for  a  Stark  or  Zeeman  shifted 
atom.  As  expected,  the  greatest  effect  occurs  with  the  linear  Stark  shift;  however, 
few  atoms  or  molecules  exhibit  such  a  shift.  The  quadratic  Stark  effect  and  the 
Zeeman  effect  have  comparable  ranges  in  which  the  Townes— Merritt  could  be 
observed.  These  are  generally  limited  by  the  maximum  field  strength  available. 
Experimental  constraints  for  level  crossings  and  quantum  beats  are  similar  to  the 
Townes-Merritt  effect  calculations. 

Conditions  for  observing  quantum  beats  in  linear  and  quadratic  Stark  shifted 
atoms  are  similar  to  those  for  the  Townes-Merritt  effect  previously  discussed.  To 
show  this,  consider  the  power  radiated  by  a  three  level  atom  in  an  RF  field.  For  the 
linear  Stark  case,  this  is  given  by  equation  (4.34) 


P  =  2P0e 


-At 


1  + 


luf)  cos(w]2-m/)t 


(5.14) 


The  signal  decays  exponentially,  at  rate  A,  with  a  superimposed  harmonic 
modulation  defined  by  the  terms  in  the  Bessel  function  series.  Oscillations  occur 
with  frequency  w,2— nv  and  relative  amplitude  Jn(^)  •  The  modulation  factor  Mn 
for  a  given  component  is  in  fact  equal  to  the  Bessel  function  coefficient. 
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For  the  beats  to  be  observable,  two  conditions  must  be  met:  the  modulation 
factor  has  to  be  large  enough  to  be  detected  and  the  beat  frequency  must  be  within 
the  measuring  capability  of  the  experiment.  Much  as  was  discussed  in  the  Townes- 
Merritt  effect,  to  detect  beats  resulting  from  interactions  with  the  sidebands,  the 
Bessel  function  argument  should  be  at  least  of  order  unity,  so  higher  order  terms  are 
significant.  The  conditions  on  £lv  and  v  are  the  same  here  as  for  the  linear  shift  in 
the  Townes— Merritt  effect,  resulting  in  a  wide  range  of  possible  transitions. 

The  second  consideration  is  the  beat  frequency,  and  the  primary  concern  is 
the  signal  bandwidth  in  the  experiment.  For  example,  in  a  decay  rate  measurement, 
the  sample  is  repeatedly  excited,  then  the  radiated  power  is  measured  during  a  finite 
time  period  At  which  is  delayed  from  the  initial  excitation.  As  the  delay  time  is 
increased,  the  power  received  decreases  according  to  the  decay  curve.  The  time 
period  At  establishes  a  maximum  bandwidth,  and  hence,  sets  a  limit  on  the 
measurable  beat  frequency. 

Given  a  limiting  bandwidth,  BW,  for  a  particular  experimental  setup,  only 
those  frequencies  which  are  less  than  that  bandwidth  are  detectable;  the  others 
average  to  zero.  For  example,  assume  BW  =  50 A  and  ui12  =  100A.  If  the  RF 
frequency  is  on  the  order  of  30A  ,  then  the  only  measurable  beats  are  those  whicn 
satisfy 

1 100A  -  n30A  |  <  50A  (5.15) 

In  this  case,  n  =  2,3,4  , corresponding  to  beat  frequencies  of  10A,  20A,  and  40A. 

It  would  be  possible  to  adjust  the  RF  frequency  and  the  bandwidth  such  that 
only  one  beat  signal  is  present,  and  the  power  radiated  would  be 

P  =  2P 0e~At  1  +  Jn(^)  cos(u/12— m/)t|  (5.16) 

j 
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This  raises  the  possibility  of  using  the  beat  signal  to  measure  the  separation  between 
widely  separated  upper  states,  with  a  large  ul2  ■ 

There  is  one  further  observation  concerning  the  beat  signals.  It  is  possible 
that  the  beat  signals  resulting  from  interaction  with  the  RF  field  could  be  obscured 
by  the  zero  field  beat  signal,  ul2.  This  zero  field  signal  can  be  eliminated  by 
adjusting  the  field  strength,  such  that  the  Bessel  function  coefficient  J0(^*')  =  0  . 
The  first  zero  occurs  for  ^  =  2.4. 

In  the  case  of  the  quadratic  Stark  shift,  the  experimental  considerations  are 
very  similar  to  the  linear  case.  From  equation  (4.48),  the  power  radiated  is 


P  =  P0e  At  1  +K2  +  2K^Jn(^)  cos(u;21+  §'+2n//)t  (5.17) 

n 

The  primary  differences  are  the  constant  of  proportionality  K  (defined  by 
|  /x10 1  =  K  |  fil0 1 )  and  the  beat  frequencies  occur  on  the  even  harmonics  of  v.  In 
addition,  the  Rabi  frequency  is  now  replaced  by  fP  ,  which  is  equal  to 

2  2 

ft1  =  <**“1  -  (5.18) 

U\U)2  V  7 


Assuming  u>i  »  ,  this  can  be  rewritten  as  ST'  =  —  (1-K2)  . 

The  modulation  factor  for  the  nth  beat  signal  is  then 


2K  T  ft2  „  \/2 
~  T+K?Jn 


(5.19) 


which  gives  the  relative  intensity  of  the  beat  signal  with  frequency  o>21+  j  +2nv 
If  K  «  .5  ,  then 

M"  =  5  Jn(^ 
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This  is  very  similar  to  the  Bessel  function  coefficient  in  the  quadratic  case  of  the 
Townes— Merritt  effect;  thus,  the  same  overall  conditions  for  the  experimental 
parameters  apply.  As  in  the  linear  case,  the  beat  frequencies  must  be  less  than  the 
bandwidth  to  enable  detection. 

In  a  level  crossing  experiment,  instead  of  determining  the  radiated  power  as  a 
function  of  time,  the  total  energy  radiated  in  a  specific  direction  is  measured.  As 
shown  in  the  last  chapter,  variations  in  the  total  energy  as  a  function  of  the  RF 
frequency  are  indications  of  level  crossings.  As  in  quantum  beats  and  the 
Townes— Merritt  effect,  these  variations  are  detectable  if  the  Bessel  function 
coefficients  for  higher  order  terms  are  significant.  Again,  this  requires  that  the 
argument  be  at  least  of  order  unity.  One  other  consideration  in  a  level  crossing 
experiment  is  the  observation  angle.  The  initial  calculations  were  performed  with 
peaks  indicating  the  level  crossings;  another  geometry  may  be  more  favorable  for 
precise  measurement  of  the  RF  frequency  at  which  the  crossings  occur.  For  the 
purpose  of  brevity,  only  the  first  case,  where  p,0  =  /i0t+  and  =  p0e„  ,  will  be 
used;  and  assume  the  observation  is  in  the  xy  plane,  so  0  =  ^  and  <p  may  vary. 

From  equation  (4.67),  the  total  energy  in  the  linear  Stark  shifted  case  is 


where  An  =  ul2— nv  .  For  <j>  =  0  ,  this  results  in  a  series  of  peaks  occurring  when 
An  =  0  ,  corresponding  to  a  level  crossing  condition.  The  relative  magnitude  of  each 
peak  is  so  the  detectability  then  depends  on  the  RF  field  strength  and 

frequency,  as  before.  In  fact,  the  same  requirements  exist  here  as  for  the  linear 
Stark  shifted  case  of  quantum  beats. 
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Figure  18.  Le'  el  Crossing  Signal  for  (p  =  */ 4 

For  experimental  purposes,  it  may  be  more  convenient  to  shift  the  detector 
to  <f>  =  j  .  In  this  case,  the  detector  is  no  longer  in  line  with  the  initial  excitation 
pulse,  reducing  potential  degradation,  and  the  level  crossing  signal  has  a  different 
shape.  For  a  single  level  crossing  event,  where  ni/  »  ul2 ,  equation  (5.22)  becomes 


(5.23) 


Figure  18  depicts  the  energy  detected  near  a  level  crossing;  the  crossing  occurs  when 

An  =  0  or  when  £  passes  through  one.  The  precise  RF  frequency  at  which  the 
J  O 

crossing  occurs  can  more  easily  be  determined  in  this  case. 

In  the  case  of  a  Zeeman  shift,  the  level  crossing  equation  takes  i  hc  same  form 
as  equation  (5.22)  except  the  Bessel  function  coefficient  is  replaced  with  Jn(^*')  . 
The  conditions  for  observing  level  crossings  here  are  the  same  as  were  discussed  in 
the  Townes— Merritt  effect. 
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For  a  quadratic  Stark  shift,  the  total  energy  radiated  is  listed  in  table  6.  For 


1 


case  I, 


}=j« 


1  +  K‘  + 


A2cos  2  <j> 

+ 

AAn'sir  2 <j> 

"" FJ 

+ 

n  J 

(5.24) 


O' 

A  level  crossing  occurs  when  =  u>i2  +  -j-  +  2rw  =  0  .  The  relative  intensity  of 
the  individual  peaks,  when  <p  =  0  ,  is 


\J  _  K2  ,ft'. 

[7o]max'T+T^Jn(^ 


If  K  a  l/|i?  ,  then  using  equation  (5.18),  this  becomes 

J  _  2  T  (  fi2 
Jo.  max  5  11  \u}0v 


(5.25)  ' 


(5.26) 


which  is  similar  to  equation  (5.20)  in  the  quadratic  case  of  quantum  beats.  Thus,  in 
order  to  detect  «.t  least  the  first  order  peak  (n=l)  ,  the  minimum  RF  field  strength 
is  about  1000  V/cm  for  an  atomic  transition  frequency  less  than  5  *  1014  Hz  and  v 
~  107  Hz. 

These  order  of  magnitude  calculations  demonstrate  the  oossibility  of 
detecting  quantum  beats  and  level  crossings  resulting  from  interaction  between  the 
atom  and  a  nonresonant  RF  field.  Thus,  it  is  feasible  to  use  the  level  crossing 
technique  to  measure  the  spacing  between  nearby  excited  states,  making  this  a 
useful  spectroscopic  tool. 
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VI. 


The  primary  purpose  of  this  research  was  to  investigate  the  possibility  of 
using  a  nonresonant  electromagnetic  field  as  a  means  to  measure  closely  spaced 
atomic  energy  levels  by  causing  a  level  crossing.  Level  crossing  experiments  have 
been  used  to  make  spectroscopic  measurements;  however,  in  these  experiments, 
static  electric  or  magnetic  fields  were  used  to  give  rise  to  the  energy  shifts.  Using 
static  fields  has  the  disadvantages  of  requiring  relatively  strong  fields  to  create  a 
significant  shift,  field  nonuniformities  result  in  broadening  and  lower  resolution,  and  . 
the  lack  of  precision  in  measuring  the  actual  field  strength  leads  to  significant 
uncertainty.  If  an  oscillating  field  is  used  to  create  the  energy  shifts  and  subsequent 
crossings,  these  problems  are  alleviated. 

A  nonresonant  field  causes  an  effective  splitting  in  the  atomic  energy  states 
in  what  is  referred  to  as  the  Townes-Merritt  effect.  The  splitting  can  be  viewed  as 
individual  states  resulting  from  the  combination  of  atom  and  field  states  into  a 
single  quantum  system.  In  this  'dressed  states'  picture,  photon  ladders  are 
superimposed  on  the  atomic  energy  states.  Transitions  between  the  new  states 
result  in  sidebands  on  the  original  atomic  line  separated  by  multiples  of  the 
nonresonant  field  frequency.  The  Townes-Merritt  effect  for  a  two  level  atom  was 
analyzed  here  with  a  quantum  electrodynamic  (QED)  approach,  using  the 
Heisenberg  operator  formalism  to  describe  the  atomic  dynamics.  The  adiabatic 
approximation  was  used  to  include  the  effects  of  the  nonresonant  field  on  the  atomic 
energy  states  for  both  linear  and  quadratic  Stark  shifts.  The  description  of  the 
Zeeman  effect  is  mathematically  equivalent  to  the  linear  case.  The  results  agree 
with  the  dressed  states  model  and  the  original  experiment. 

In  a  multi-level  atom,  the  Townes-Merritt  splitting  can  result  in  other 
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effects:  specifically  quantum  beats  and  level  crossings.  These  effects  result  from 
interference  between  separate  transitions  which  originate  in  nearby  excited  states 
and  end  in  a  common  ground  state.  Level  crossings  ana  quantum  beats  were 
analyzed  for  linear  and  quadratic  Stark  shifted  atoms,  again  using  the  Heisenberg 
operator  formalism.  Quantum  beats  which  depend  on  the  field  frequency  were 
predicted  in  both  cases,  with  the  stronger  effect  in  the  linear  case.  If  the  quantum 
beat  signals  can  be  detected,  this  may  provide  a  means  to  determine  accurately  the 
difference  between  widely  spaced  energy  states. 

In  the  case  of  level  crossings,  expressions  were  developed  defining  the  total 
radiated  energy,  along  a  specific  viewing  axis,  as  a  function  of  the  field  frequency.  A 
level  crossing  results  in  a  change  in  the  fluorescence  pattern  indicated  by  a  variation 
in  the  received  energy.  The  peaks  in  the  energy  signal  are  relatively  narrow,  on  the 
order  of  2A,  allowing  precise  measurement  of  the  crocsing  point.  This  can  then  be 
used  to  accurately  determine  the  original  energy  spacing.  Since  the  energy  shifts 
here  are  a  function  of  the  field  frequency,  using  an  oscillating  field  alleviates  the 
requirement  for  a  strong  field  in  the  case  of  a  static  shift.  In  addition,  a  nonuniform 
field  does  not  result  in  broadening  as  in  the  static  case,  and  the  critical  parameter, 
frequency,  can  be  measured  with  extreme  accuracy.  Thus,  using  a  nonresonant  field 
in  a  level  crossing  experiment  can  result  in  a  very  accurate  spectroscopic  technique. 

Experimental  confirmation  of  this  technique  will  be  the  final  proof  of  its 
usefulness.  The  last  chapter  provided  a  look  at  some  of  the  critical  parameter  in 
such  an  experiment.  It  was  shown  that  for  an  idealized  atom,  there  is  a  wide  range 
of  possible  transition  frequencies  that  are  applicable,  infrared  to  the  mid— visible. 
Typical  RF  field  frequencies  are  on  the  order  of  107  to  1010  Hz  and  field  strengths  of 
100  to  1000  V/cm.  The  obvious  follow-on  work  is  to  find  an  appropriate  atomic  or 
molecular  sample  and  conduct  an  experiment  to  confirm  this  technique. 
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